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AGE-0 SPOT (LEIOSTOMUS XANTHURUS) FROM TWO 
ESTUARIES ALONG CENTRAL FLORIDA'S EAST 
COAST: COMPARISONS OF THE TIMING OF 
RECRUITMENT, SEASONAL CHANGES IN 
ABUNDANCE, AND RATES OF GROWTH 
AND MORTALITY. 


RICHARD PAPERNO(“)? 


Florida Fish & Wildlife Conservation Commission, Florida Marine Research Institute, 
Indian River Field Laboratory, 1220 Prospect Ave, Suite 285, Melbourne, Florida 32901 


ABSTRACT: Data from fixed-station sampling conducted during 1991-1997 were used to 
determine the timing of recruitment, the seasonal changes in abundance, and the growth and 
mortality rates of age-O spot (Leiostomus xanthurus) populations in the Indian River Lagoon 
(IRL) and Halifax River (HR) along the Atlantic coast of Florida. The magnitude of recruit- 
ment was greater in the HR system than in the IRL, particularly during the comparatively 
weak year-class of 1994. In the IRL, juvenile growth rates varied between years, with the 
instantaneous growth coefficient (G) ranging from 0.28 to 0.73. In the HR, however, growth 
varied little, with G ranging from 0.36 to 0.43. Estimates of mortality (Z) ranged from 0.02 
to 0.08/d for the IRL and from 0.02 to 0.03/d for the HR. For age-O spot, recruitment, growth, 
and mortality varied greatly between years and systems and may be strongly constrained by 
the impact advective processes have on larval supply to the estuary and by the distance to 
the nursery habitats has on the distribution and survival of settling juveniles within the 
estuary. 
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SPOT, (Leiostomus xanthurus Lacépede), is an often-studied member of 
the economically important family Sciaenidae. The species ranges from Mas- 
sachusetts on the Atlantic coast of the United States to Mexico in the Gulf 
of Mexico and is a seasonally abundant inhabitant of estuarine waters (Hales 
and Van Den Avyle, 1982). Spot have a protracted spawning season that 
lasts from October through May. Peak spawning occurs offshore in winter 
(Dec. and Jan.), and juveniles subsequently recruit to shallow inshore areas 
(e.g., tidal creeks and seagrass flats). The status of the commercial harvest 
has been studied in Florida’s coastal waters (Muller et al., 1996; McRae et 
al., 1997), but little is known about the early life history of spot, and no 
attempt has been made to estimate growth and mortality rates of age-O fish 
during their residency in Florida estuaries. 

Many aspects of the life history of spot from estuaries in the mid-At- 
lantic and the northern Gulf of Mexico have been studied (e.g. Sundararaj, 
1960; McCambridge and Alden, 1984; O’Neil and Weinstein, 1987), but 
little attention has been given to populations along Florida’s Atlantic or gulf 
coasts (Townsend, 1956; Springer and Woodburn, 1960). In the present 
study, I compare the timing of recruitment, the seasonal changes in abun- 
dance, and annual estimates of growth and mortality rates of two estuarine 
populations of age-O spot along the Atlantic coast of Florida, USA. 


MATERIALS AND METHODS—Spot were collected from two estuaries along Florida’s Atlantic 
coast: the Halifax River estuary (HR) and the Indian River Lagoon (IRL; Fig. 1). The HR 
sampling area consisted of the Halifax River from Ponce de Leon Inlet north to the Tomoka 
River (29° 24’ N) and south to include the Mosquito Lagoon (28° 44’ N; Fig. 1A). The IRL 
sampling area consisted of the northern Indian River Lagoon (north of 28° 01’ N), excluding 
the Mosquito Lagoon subbasin (Fig. 1B). In both systems, sampling stations were widely dis- 
tributed geographically and were located in several different habitat types. In the HR, 27 stations 
(13 seine, 14 trawl) were sampled monthly from January 1993 to September 1997 (56 months). 
During January 1993 to December 1995, five consecutive seine or trawl hauls were made at 
stations | through 16, and three consecutive hauls were made at stations 17 through 27. In 
1996 and 1997, three consecutive hauls were made at each of the 27 stations. Stations 1 through 
8 were located approximately 30 kilometers to the north of Ponce de Leon Inlet in and around 
the Tomoka River. This area has a mud and detritus bottom, no submerged aquatic vegetation, 
and is influenced by tides (mean range: 0.2 m) and by freshwater input from the Tomoka River. 
Stations 9 through 16 were located in the vicinity of Ponce de Leon Inlet (Fig. 1A). The 
estuarine area around Ponce de Leon Inlet region is characterized by no submerged aquatic 
vegetation and a sand bottom that is heavily influenced by inlet dynamics (mean tidal range of 
0.7 m; Taylor and Yanez, 1992). Stations 17 through 27 were located to the south of Ponce de 
Leon Inlet in the Mosquito Lagoon and ranged between 15 and 45 kilometers from the inlet. 
Mosquito Lagoon is a relatively shallow subbasin of the Indian River Lagoon with a vegetated 
mud-sand bottom. Tidal influence is negligible, and water circulation is predominantly wind 
driven. In the IRL, 20 stations (14 seine, 6 trawl) were sampled monthly from April 1991 to 
March 1996 (62 months). Stations in the IRL were located north of Sebastian Inlet and ranged 
from approximately 15 to 80 kilometers from the inlet. Both systems contain salt-marsh and 
spoil-island habitats that have shore vegetation comprised of salt marsh cordgrasses (Spartina 


' richard.paperno @ fwe-.state.fl.us 


No. 2 2002] PAPERNO—EAST COAST AGE-0 SPOT 87 


Atlantic 


Ocean 
_Helifax River 


Atlantic 


Ocean 


Cape Canaveral 
Kilometers 
i ecalessec sl 


0 25 50 


indian River 


Kilometers 
T 
100 


Fic. 1. Map of Florida showing the locations of the Halifax River-Mosquito Lagoon and 
Indian River Lagoon. 


spp.), glasswort (Salicornia spp.), saltwort (Batis maritima), and mangroves (Rhizophora man- 
gle, Avicennia germinans, Laguncularia racemosa). 

In both estuaries, a 21.3-m center-bag seine with 3.2-mm mesh was deployed along shore- 
lines in water depths less than 1.5 m. Gently sloping shorelines were sampled by manually 
dragging the seine 9.1-m perpendicular to the shore and then hauling the net onshore. Shallow 
offshore seagrass flats were sampled by manually dragging the seine 9.1-m and then hauling 
the poles together and retrieving the net offshore. Steeply banked shorelines were sampled by 
deploying the seine along the shore in an arc from the stern of the boat. In waters deeper than 
1.8 m, a 6.1-m otter trawl with 3.2-mm mesh liner in the codend was towed at approximately 
3.7 km hr~' for five minutes (avg. distance = 0.185 km/sample). 

For all hauls, fish were identified to the lowest taxa, enumerated, and up to 100 fish/taxa 
were measured to the nearest mm standard length (SL) prior to release. Large samples (>1000 
individuals) were subsampled with a modified Motoda box splitter (Winner and McMichael, 
1997), and the total number was estimated from the subsample. At each station, water temper- 
ature [°C] and salinity [%0] measurements were collected with a Hydrolab®Surveyor II multi- 
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probe. Water depth, bottom type, and a visual assessment of the percentages of seagrass cover 
(for seine sites) were recorded at each site. All data were collected during daylight hours. 

Monthly abundance of age-0 spot was estimated at the seine sites in each system; these 
estimates of abundance were used to describe annual variations in the timing and magnitude 
of recruitment and seasonal changes in abundance of age-O spot in each of the estuarine systems. 
The mean numbers of age-O spot per 100 m?* were calculated. To ensure that only age-O spot 
were included in the various analyses (abundance, growth and mortality) and to reduce bias 
resulting from incomplete recruitment to the gear and emigration of larger fish, the data included 
only spot between 20 to 60 mm SL collected during the recruitment window of January to 
April. The size limit and recruitment window were selected to account for the continuous 
recruitment of new smaller individuals and loss of larger emigrating individuals by comparing 
the length-frequency distributions of spot at the shallow- and deep-water sites in the two es- 
tuaries (Kolmogorov-Smirnov two-sample test; Sokal and Rohlf, 1981). Monthly length-fre- 
quencies based on proportions of fish collected in each length class were combined over years 
to describe within-year variability rather than between year variability within each system. 
Pearson product moment correlations were used to determine whether seasonal variations in 
age-O spot abundances were correlated with mean surface temperatures. The annual variations 
in age-O spot abundance were examined at seine sites by year, bottom type (mud or sand), 
percent cover of seagrass, month, and salinity, and the temperature effects by estuary. General 
linear models and PROC GLM (SAS Institute, 1988) were applied to In(x+1) transformed data 
prior to analysis to generate standardized monthly and annual abundance indices from the first 
haul at each site. Comparisons of the annual abundance indices were obtained with Tukey’s 
multiple comparison test (Sokal and Rohlf, 1981). 

Length-frequency data from seine collections were used to calculate growth rate estimates 
for each year class of spot in each estuary. Spot initially recruit to shallow-water areas and 
remain in these areas through the early juvenile phase of their life history (Weinstein and 
Walters, 1981). Therefore, to minimize bias associated with the emigration of larger age-O spot, 
only data from seine collection were used. Without the adjustment for emigration, the growth 
estimates would tend to be biased toward the slower-growing (smaller) individuals within the 
population, and therefore the mortality estimates would be inflated. Instantaneous growth co- 
efficients (G) were calculated by using the average time that it took fish to grow to a certain 
length in each estuary during each year (DeAngelis et al., 1980). This method utilizes a partial 
differential equation model to account for the nonuniformity of individual lengths that result 
from protracted recruitment of individuals over a long period of time. Growth rates were derived 
(Eqns): 


InL, = InLy + G X (#/30), (1) 
where G = instantaneous growth coefficient; L, = monthly mean length (mm SL); Lp = theo- 
retical size at which spot recruit to each estuary; and (t/30) = time in days. 

Mortality rates were estimated by combining length-frequency and abundance data from 


seine collections made in each estuary during each year. Mortality rates were estimated (Eqn. 
2): 


In (n,) = a +m X In (J), (2) 


where n, = the total number of spot at length 1; / = length class 1; m = slope. The instantaneous 
daily mortality (Z) could then be estimated (Eqn. 3): 


Z = {G(m + 1)}/30, (3) 


where Z could be estimated from G and m (Deegan, 1990). To reduce the bias that results from 
incomplete recruitment of smaller individuals to the gear or emigration of spot as they grow 
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Fic. 2. Monthly indices of relative abundance for age-O0 spot and mean temperature from 
the Halifax River (1993-1997) and Indian River Lagoon (1991-1996). 


(as described above for growth rate estimates), only fish between 20 and 60 mm SL were used 
for mortality calculations. 


RESULTS—Seasonal patterns of the timing of recruitment were found to 
be similar in the IRL and the HR systems (Figs. 2—4). In each system, spot 
began to be fully recruited between 15 to 20 mm SL in January as temper- 
atures approached seasonal lows (18—20°C in the IRL and 15—18°C in the 
HR; Figs. 2-4). Annually, 1993 represented strong recruitment years in both 
systems. In addition, the magnitude of recruitment in the respective systems 
was different during the strong 1993 and weak 1994 year classes, with more 
spot being collected on average from the HR system than from the IRL (Fig. 
Dy: 

In the HR system, similar-sized age-O spot were collected in seines and 
trawls from January to April (K-S test, P>0.05; Fig. 3). However by April, 
although few in number, spot >60 mm began to appear in trawl samples 
from the IRL (K-S test, P = 0.05), but not in those from the HR (K-S test, 
P>0.05). After April, spot larger than 60 mm were collected in trawls in 
both systems, although in the IRL, spot of any size were completely absent 
or rare in samples collected after July (Figs. 3 and 4). In HR from May 
through December, the proportion of age-O spot >60 mm collected in trawls 
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Fic. 3. Length-frequency distributions of spot captured at seine ( ) and trawl (-:------ ) 
stations in the Halifax River-Mosquito Lagoon from 1993 to 1997. N, is the number of fish 
measured from seine samples. N; is the number of fish measured from trawl samples. 


was greater than those from seines (Fig. 3). Although the processes of im- 
migration and emigration in the shallow-water areas (seine stations) can bias 
the rate if decline in abundance was used to estimate mortality, only data 
from months over which these processes appeared low (January—April) were 
used for subsequent growth and mortality estimates. 

Juvenile spot abundance varied significantly (P<0.05) between years in 
both estuarine systems. The final GLM models each had year as the only 
main effect with 7? = 0.088—0.139 (Table 1). Annual spot abundance was 
generally greater in HR (2.4—28.2 fish/haul) than in the IRL (0.6—20.8 fish/ 
haul; Fig. 5). There was also a difference in the variability associated with 
the annual abundance estimates. The 95% confidence intervals around the 
annual abundance estimates were smaller for the HR (< | fish per set) than 
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Fic. 4. Length-frequency distributions of spot captured at seine ( ) and trawl (--::---: ) 
stations in the Indian River Lagoon from 1991 to 1996. N, is the number of fish measured 
from seine samples. N; is the number of fish measured from trawl samples. 


those for the IRL (Fig. 5). In both systems, 1993 and 1995 were strong 
recruitment years, whereas, 1994 was a poor recruitment year for juvenile 
spot (Fig. 5). Spot abundance in the HR from 1994 through 1997 was rel- 
atively stable but significantly lower than during 1993, which was a strong 
recruitment year. Abundance of age-O spot in the IRL had a significant neg- 
ative correlation (r = —0.26, P<0.05) with temperature; fish were essentially 
absent during the months of September to December (Figs. 2 and 3). In the 
HR, the relationship was also negatively correlated (r = —O.21) but not as 
Statistically significant (P = 0.11), and fish were collected more consistently 
throughout the year (Figs. 2 and 4). 

Year-class growth rates varied between years and followed no discern- 
ible pattern in the two estuaries. In the HR, the instantaneous growth co- 
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Fic. 5. Indices of relative abundance for age-O spot in the Halifax River (1993-1997) 
and Indian River Lagoon (1991-1996). Error bars represent 95% confidence intervals. Values 
with the same letter are not significantly different at P<0.05. 


efficient varied from 0.36 to 0.43/month (Table 2). In the IRL, growth rates 
varied greatly, with the instantaneous growth coefficient ranging from 0.28 
to 0.73/month (Table 2). 

Mortality rates (Z) were estimated for each year for each system. In each 
system the greatest mortality occurred during the strong recruitment year of 
1993 (Table 3). As with growth, Z varied little between years in the HR 
(O.02—0.03/d) but varied greatly in the IRL system (0.02—0.08/d). 


DIscUssION—Spot, which have a spawning season that can last from 
October through May, recruit to shallow estuarine habitats during winter and 
early spring and may reside in the estuaries until fall. Although the onset of 
recruitment begins earlier in Florida than in more temperate estuaries 
(Springer and Woodburn, 1960), the general seasonal pattern of age-O spot 
occurrence in the IRL and HR was consistent with the seasonal patterns of 
occurrence reported for spot in other systems along the western Atlantic and 
Gulf of Mexico (Sundararaj, 1960; Chao and Musick, 1977; McCambridge 
and Alden, 1984). Although the timing of recruitment events was found to 
be similar in both systems, the magnitude of the recruitment events was 
generally greater in the HR system. This may be due in part to the location 
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TABLE 2. Regression statistics for mean standard length (SL) versus month (m) for age-0 
spot. The regression equation takes the form: Ln (SL) = In(Z,) + G X m where In(Z)) is the 
intercept and G is the instantaneous growth rate (slope). Only data for fish less than 60 mm 
SL were used. All slopes and intercepts were significantly different from zero (P < 0.05). 


Year In(Lo) SE G SE li 
Indian River Lagoon 
1992 1.249 0.111 0.730 0.055 0.956 
1993 1.982 0.182 0.435 0.062 0.874 
1994 2.669 0.116 0.276 0.039 0.891 
1995 Dae 0.090 0.394 0.028 0.966 
1996 opie 0.107 0.408 0.030 0:95) 
Halifax River Estuary 
1993 2.156 0.044 0.363 0.014 0.989 
1994 1927) 0.112 0.434 0.045 0.931 
1995 21037; 0.040 0.420 0.015 0.991 
1996 Qalis 0.102 0.360 0.033 0.943 
1997 1.801 0.109 0.432 0.043 0.935 


of Ponce de Leon Inlet within the HR sampling area, to differences in larval 
predation pressure between the areas, to differences in larval supply, or 
merely to a characteristic of the sampling sites. These four possibilities, 
singly or in combination, could result in differences in local population 
densities of settling juveniles. 

Differences in larval supply of spot to the two systems may be due to 
a difference in the magnitude of spawning events or supply of spot that 
support each system or to less favorable conditions (wind, current, distance 


TABLE 3. Regression statistics for total number versus length (SL) for age-O spot. The 
regression equation takes the form: Ln (total number) = a + m x In(length) where a is the 
intercept and m is the slope. Only data for fish less than 60 mm SL were used. Estimates of Z 
are calculated with equations in text. 


In(total 

Year number) SE a SE r Jb, 
Indian River Lagoon 

1992 12.742 ah63 — 1.647 0.318 0.817 0.0644 

1993 12.907 0.482 —4.585 0.368 0.963 0.0810 

1994 7.260 5.645 — 6.862 1.426 0.268 0.0171 

1995 8.104 1.506 AU) 0.412 0.394 0.0238 

1996 12.599 0.696 =202i 0.227 0.930 0.0411 
Halifax River Estuary 

1993 12.656 0.763 — 1.654 0.209 0.913 0.0322 

1994 8.946 0.682 =O09873 0.187 0.785 0.0271 

1985 8.581 0.718 —0.662 0.193 0.703 0.0232 

1996 8.740 0.468 =O Ae 0.128 0.858 0.0212 


ey 5.852 1.485 —0.256 0.406 0.063 0.0181 
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to nursery area) for larval recruitment. Significant relationships have been 
demonstrated between larval supply and onshore and longshore winds, and 
current velocity and direction that suggests the magnitude and distribution 
of the larval supply to and within a nursery area may be related to complex 
interactions with the physical environment (Thorrold et al., 1994; Tilney et 
al., 1996; Epifanio and Garvine, 2001). Although Ponce de Leon Inlet (HR) 
has a more energetic tidal exchange than the Sebastian Inlet (IRL), condi- 
tions at Sebastian Inlet (net inflow) and in the intracoastal within the north- 
ern Indian River Lagoon appear to support transport to the northern reaches 
of the lagoon during late fall and early spring (Taylor and Yanez, 1992; 
Smith, 2001). It is also possible that internal northern transport mechanisms 
may not be sufficient in some years to overcome the difference in distance 
between the respective inlets and nursery grounds. Sebastian Inlet, which 
serves as the connection to the Atlantic for the northern IRL, is located 
approximately 15 km from the southern extremity of the IRL study area, 
whereas Ponce de Leon Inlet is located near the center of the HR study area. 
No evidence has been presented that indicate that spot spawn in estuaries, 
so this difference in distance to the respective nursery may have a significant 
effect on the comparative magnitude and variability of annual recruitment 
for spot and other estuarine-dependent species that spawn offshore, near HR 
and IRL. Since conditions in the northern IRL should support immigration 
into and distribution within the system, differences in recruitment between 
the HR and IRL systems may be a result of differences in advective transport 
of larvae from nearshore spawning areas to the estuaries and for spot may 
be compounded by distance to the nursery grounds. Differences in advective 
transport have been suggested to be the cause of variation in recruitment of 
bluefish to estuaries in the Mid-Atlantic Bight (McBride and Conover, 1991; 
Epifanio and Garvine, 2001). 

Although the general seasonal timing of juvenile spot recruitment in 
Florida was found to coincide with the timing of spot recruitment in other 
areas, spot in Florida differed in cohort development. Evidence of distinct, 
multiple spawns was found in North Carolina (Weinstein and Walters, 1981), 
South Carolina (Beckman and Dean, 1984), and Virginia estuaries (Mc- 
Cambridge and Alden, 1984) but not in Gulf of Mexico or Florida coastal 
systems (Sundararaj, 1960; this study). In Florida estuaries, juvenile spot 
recruitment is a continuous process that can last from January through April. 
In general, the presence of multiple cohorts of juveniles throughout tem- 
perate nursery habitats may enhance year-class strength in the face of en- 
vironmental conditions that occasionally extend beyond the range of tolerance 
of a particular cohort (Lambert, 1984; Lambert and Ware, 1984; Szedlmayer 
et al., 1990; Paperno et al., 2000). In Florida, however, environmental con- 
ditions (e.g., warm temperate climate) may not be such that this strategy is 
necessary to maintain the integrity of the year class, and other factors, such 
as predation rate or advective processes, may be more important in deter- 
mining juvenile year-class strength (Epifanio and Garvine, 2001). 
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Distribution within and recruitment of fish to an estuary may be influ- 
enced by factors such as water temperature, salinity, food availability, pre- 
dation, competition, and habitat structure or availability (Sogard, 1992; Milt- 
ner et al., 1995; Jones et al., 1996; Young et al., 1997). Age-O spot (<60 
mm SL) are ubiquitous within the IRL and HR systems even during poor 
years of recruitment. Juveniles inhabit both shallow-water habitats (<1.8 m) 
and deeper channels (>1.8 m), although they are more abundant in the 
shallows. Movement of fish out of the shallow-water areas to the channels 
was evident earlier in the IRL than in the HR system. The HR is a deeper 
water system than the IR is and has a smaller amount of shallow-water 
seagrass and sand-flat habitat. In the HR system, the deep water near shore 
may have a moderating effect on the water temperature there, which may 
explain why spot in the HR system move to deeper water later than spot in 
the IRL system do. In a thermally affected estuary, temperature has been 
demonstrated to have a direct effect on the distribution of some estuarine 
fishes: fish densities can decrease during warmer months and increase during 
cooler months (Jones et al., 1996). Thus lower (~4°C) temperatures in late 
winter in the HR system, along with the potential deep water refugia, may 
be the principal reason that spot there remain in the estuary later in the 
season than spot in the IRL do. 

Year-class growth rates of age-O spot varied between years and followed 
no discernible pattern in the two estuaries. Growth rates estimated for spot 
from both systems were comparable with rates of spot in other estuarine 
systems. To compare growth rates for age-O spot from this study (G = 0.28— 
0.73) with those found from a more northern estuary, the same growth equa- 
tion was applied to spot data presented by Weinstein and Walters (1981) for 
spot collected from the Cape Fear estuary, North Carolina. Instantaneous 
growth rates for spot from Cape Fear estuary (G = 0.23—0.31) overlapped 
at the lower range with those from this study. This is not surprising when 
one considers that in 1994, the year that the lowest growth rates were found, 
was also the year of the lowest monthly mean temperatures. Thus, conditions 
in the northern Indian River Lagoon may have been similar to those during 
the Weinstein and Walters (1981) study. 

Age-O spot experienced similar rates of mortality in both east-central 
Florida systems during most years, but as with the growth rates, the inter- 
annual variation in the mortality rates was less in the HR than in the IRL. 
This difference in interannual stability may be a reflection of the relative 
influence of local environmental conditions on survival of juvenile fishes 
residing in a tidally driven system (HR) versus a shallow, wind-driven sys- 
tem (IRL). Daily mortality rates calculated for spot in this study (Z = 
—0.017 to —0.081) overlapped in range with mortality rates calculated for 
spot in several mid-Atlantic estuaries (Z = —0.023 to —0.041) (Weinstein 
and Walters, 1981; Weinstein, 1983; Currin et al., 1984; Miller et al., 1985). 

The patterns observed in the rates of growth and mortality in age-O spot 
from these two estuarine systems may indicate that the tradeoffs between 
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growth and survival necessary for success of a species may be different for 
fish in Florida estuarine systems than for fish in more temperate estuaries. 
For many species in temperate systems, optimal growth rates and minimal 
mortality rates are maintained through habitat selection or timing of spawn- 
ing that result in reduced overwintering mortality (Sogard, 1997). Weinstein 
and Walters (1981) found that when temperatures were higher, spot ate less 
and thus grew at a slower rate than when temperatures were cool. As sug- 
gested for other species, the slower growth rates of spot in the southern 
estuaries may reflect the greater physiological constraints imposed by higher 
temperatures on southern populations (Conover and Present, 1990). In ad- 
dition, the fact that northern and southern populations have similar mortality 
rates while they reside in nursery habitat may also suggest southern popu- 
lations do not need to attain greater overwintering size because winter con- 
ditions are considerably less harsh in south Atlantic estuaries than in north 
Atlantic estuaries. 
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RE-EVALUATION OF OFF-ROAD VEHICLE 
MANAGEMENT IN BIG CYPRESS NATIONAL 
PRESERVE, FLORIDA, AFTER TWENTY-SIX 
YEARS OF MANAGEMENT 
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ABSTRACT: Jn 2000, Big Cypress National Preserve, Florida, took a large step toward 
managing off-road vehicles (ORVs) for the preservation of natural resources, culminating in 
publication of a new ORV management plan. This paper describes events and conditions that 
led the National Park Service to develop the new plan. Photographic images, trail mapping, 
hydrologic conditions, legislation, and the advent of new technologies are discussed as ratio- 
nale for restricting the level that ORVs have traditionally been allowed in the first twenty-six 
years of the preserve’s existence. 


Key Words: Big Cypress National Preserve, off-road vehicles, ecosystem 
protection 


THERE has been growing awareness regarding the suitability of off-road 
vehicle (ORV) use within units of the National Park system—including 
snowmobiling in Yellowstone National Park (Wyoming, Montana, Idaho), 
dune buggies in Mohave National Preserve (California), and most recently, 
swamp buggies and airboats in Big Cypress National Preserve (Florida). 

Big Cypress National Preserve (295,245 ha) comprises the eastern third 
of the Big Cypress Swamp in south Florida, USA. It was established in 1974 
for the purpose of protecting the upstream watershed that is vital to western 
Everglades National Park and protecting the fragile Big Cypress Swamp 
from development threats (Fig. 1). The preserve is extremely flat (slope = 
3 cm/km) and covered at the surface with a thin layer of calcium-rich soil 
called marl underlain by a limestone caprock (O—1 m depth). On a local 
scale, topographic variation as small as 60 cm separates uplands from wet- 
land (Duever et al., 1979). This local-scale variation also tends to match the 
distribution of major vegetation types in the preserve—including upland 
pinelands and hardwood hammock, herbaceous prairie, cypress strands and 
domes, and marshes. An average of 133 cm of rain falls per year, sixty 
percent of which occurs from June to October (Klein et al., 1970). Poor 
drainage results in shallow inundation of the landscape during this period. 
Water levels gradually recede during the non-rainy season (November— 
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April). All but the deepest water holes become dry by mid-winter or early 
spring, depending on the extent of rain received during the previous rainy 
season. 

Types of ORVs used in the preserve include swamp buggies, airboats, 
street-legal four-wheel drive vehicles, and all-terrain vehicles (ATVs). 
Swamp buggies have oversized tires, heavy-duty construction, and high- 
torque engines that provide adequate traction across mudded and flooded 
areas (Fig. 2). Airboats are smooth, flat-bottomed boats that operate in ex- 
tremely shallow water by means of a large propeller (2—5 meter width) that 
is fixed and caged on the back of the boat (Fig. 3). Typically, the propeller 
and engine that drive it are taken from airplanes, thus giving the boats their 
name. Street-legal four-wheel drive vehicles are also used along elevated- 
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Fic. 2. Photograph of a typical swamp buggy. Dimensions are approximately 3 m wide, 
3-5 m long, and 3-5 m tall. 


grade roads, but are not as common. Manufacturer-produced all-terrain ve- 
hicles (ATVs) have become more popular in recent times. 


DISCUSSION—Origination of ORVs in the Big Cypress Swamp—The Big 
Cypress Swamp has historically been a remote and roadless environment. 
Use of ORVs in the Big Cypress predated establishment of the preserve. 
Their use first began in the 1920s in conjunction with construction of the 
primary roads and railroad grades that bisect the interior swampland of the 
Big Cypress Swamp (Tamiami Trail, State Road 29). Logging and oil/gas 
development operations became established along the roads and utilized 
ORVs to harvest natural and geologic resources. Recreational ORV use did 
not become common until the 1940s (Tebeau, 1966). 


Allowance and use pattern of ORVs in the preserve—ORV operation is 
allowed within the preserve by virtue of the Enabling Legislation (Public 
Law 93-400) through which the United States Congress established the 
preserve in 1974, but is not a mandated activity nor are specific rules put 
forward on how to manage them. The Enabling Legislation states that “. . . 
the Secretary shall develop ... such rules and regulations as he deems nec- 
essary and appropriate to limit or control the use of Federal lands and waters 
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Fic. 3. Photograph of a typical airboat. Dimensions are approximately 2 m wide, 3—5 m 
long, and 3—5 m tall. 


with respect to: (1) motorized vehicles, ...’’, but does not provide further 
direction (U.S. Congress, 1973). 

An unrestricted use pattern best describes the manner in which ORVs 
operated for the first twenty-six years of the preserve. Permitted ORV op- 
erators were given discretion to choose the access points and backcountry 
trails they wanted to use. Exceptions to this included large-scale spatial 
closures that were instituted in 1974 (Loop Unit) and 1988 (Deep Lake 
Unit), establishment of a 2,500 annual permit quota in 1991, designation of 
a trail network in the Bear Island unit, and periodic temporal closures en- 
acted in response to extreme environmental conditions (Fig. 4). Water-filled 
borrow canals that lie adjacent to most elevated roadbeds also limited the 
total number of ORV access points to periodic structural or earthen bridges. 
In total, ORV operators have traditionally been able to gain access into the 
backcountry at approximately seventy locations. 

ORV use is least apparent during summer months due to extreme heat, 
high water, and mosquitoes. Peak ORV use occurs in tandem with hunting 
season (November—December), a time period that coincides with the end of 
the rainy season and end of the growing season. The end of the rainy season 
marks a time when a thin layer of surface water still covers most of the 
preserve, but is receding due to lack of rain. The preserve’s calcium-rich 
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marl soils are soft and easily disturbed when saturated, and provide poor 
traction to wheeled ORVs. Heavily used trails often become excessively 
mudded and rutted out, thereby making them difficult to traverse. ORV 
operators tend to avoid these areas in favor of more vegetated areas, typically 
at the sides of existing trails (causing trails to widen) or by developing new 
trails. The end of the growing season marks a time when groundcover veg- 
etation has ceased growing for the year, but is also at its maximum density. 
Duever and co-workers (1979) noted that occurrence of ORV activity at this 
time of year tends to minimize the level of ORV disturbance because veg- 
etation cover and rootmass provide additional structure to the soil, but also 
creates a condition where visual impacts caused by ORVs persist until the 
onset of the next growing season. Fig. 5 summarizes the seasonal relation- 
ship between ORV use, water levels, and growing season. 


Environmental effects of ORV use—Soils and herbaceous plant com- 
munities are most directly disturbed by ORVs use. Heavily used areas are 
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Fic. 5. Seasonal relationship among off-road vehicle (ORV) use, water levels, and grow- 
ing season in the preserve. Note that ORV use is estimated from observation since quantitative 
data is not available. 


damaged in ways that are visually apparent—including an irregular topo- 
graphic surface, excessively mudded soils, and decreased to nonexistent 
plant cover. Long-term effects in these areas are thought to include loss of 
herbaceous prairie habitat (plant loss and marl disturbance), alteration of 
surface water regimes through channelization and increased exposure of sur- 
face waters to evapotranspiration, and decreased primary production due to 
destruction of floating algal mats called periphyton. 

There has been an on-going debate regarding the extent that visual dis- 
turbance translates into long-term ecological effects. Schemnitz and Shor- 
temeyer (1974) studied the effects of airboats and track vehicles in the Ev- 
erglades, but this work was conducted on soil types not common in the Big 
Cypress Swamp. Duever and co-workers (1979) completed a resource-in- 
ventory for the preserve that described the history and pattern of ORV use 
in the Preserve. They showed an accretion of total ORV trail length over a 
three-decade period using high-altitude photography taken in 1947, 1953, 
and 1973. Stubbs (1979) described social aspects of ORV use in the preserve 
using a survey. Duever and co-workers (1981) used test plots to identify 
short-term impacts caused by various ORV and soil types. The study showed 
airboats to be the least damaging of all ORVs, and ATVs to the least dam- 
aging wheeled-ORV tested. Dwarf cypress and marl prairie were identified 
as being the most susceptible ORV-induced impacts because of their soil 
composition and relative lack of plant and root mass. Duever and co-workers 
(1986) revisited the test plots established in the 1981 study to provide an 
assessment of long-term effects. They concluded that once soils are dis- 
placed, there are few natural mechanisms capable of restoring the natural 
ground contour, and the tire ruts remain indefinitely. Welch and Madden 
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(1998) completed a database of ORV trails and high use areas. Farrell and 
co-workers (1998) described an ORV-use survey provided to off-road ve- 
hicle permit holders in the Preserve. 


Reconsideration of how ORVs are managed in the preserve—With some 
modification, the first twenty-six years of the preserve saw a continuation 
of ORV activity largely as it had existed prior to its establishment. Possible 
causes of this continued use pattern may include the vague nature of guid- 
ance provided in the Enabling Legislation, lack of information suggesting 
otherwise, an ethos that ORV-induced disturbance could be minimized by 
encouraging a dispersed-use pattern, a desire to accommodate a major user 
group of the preserve, and a tendency to keep the status quo. 

The presence of ORVs in the preserve was brought under increased 
scrutiny following a 1995 lawsuit filed against the National Park Service by 
the Florida Biodiversity Project, an environmental advocacy group. The law- 
suit contended that the level of disturbance caused by ORV use was exces- 
sive and in violation of federal statutes. A setthkement agreement between the 
two parties was reached in 1997, on the condition that National Park Service 
develop an off-road vehicle management plan that provides for improved 
management of the use. 

As part of the process to reconsider its management of ORVs, the Na- 
tional Park Service identified five factors that provide rationale for restricting 
the level that ORVs have traditionally been allowed in the preserve during 
its first twenty-six years of existence. These factors are believed to show 
that existing levels of disturbance caused by ORVs are unacceptable, will 
increase Over time if no action was taken, and contradicted legislative di- 
rection regarding ORV use provided during the establishment of the pre- 
serve. 


Photography of excessive disturbance—Aerial and ground-level photog- 
raphy highlighted unacceptable levels of soil disturbance, vegetation loss, 
and alteration of surface water inundation and flow (Fig. 6). Disturbance 
was most apparent in herbaceous prairie habitats where lack of forested 
environment allowed for unconstrained expansion and widening of trail net- 
works due to the absence of physical barriers (trees) and at major access 
points where concentrated ORV use occurred. 


Evidence of accumulated impacts over time—Mapping work performed 
by previous investigators demonstrated that the magnitude and spatial extent 
of disturbed soils has increased over time. Duever and co-workers (1979) 
reported an incremental increase in ORV trails in the 30 years prior to the 
preserve’s establishment by using high-altitude aerial photographs available 
from the USDA Soil Conservation Service for 1940 and 1953 (scale = 1: 
63,000) and Mike Hurd for 1973 (scale = 1:75,000) to map trails that were 
visible without magnification. This work estimated total trail lengths of 160, 
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Fic. 6. (a) Aerial and (b) ground level photographs of off-road vehicle damage in the 
preserve. 


250, and 1,100 km in 1940, 1953, and 1973, respectively, with the increase 
from 1953 to 1973 being explained almost exclusively as a result of recre- 
ational ORV use. Similar mapping techniques applied by NPS staff to L: 
63,000 scale aerial photography taken in 1984 estimated 1,995 km of ORV 
trails. Duever and co-workers (1979) noted that additional trails would be 
visible through magnification or use of aerial photography taken at lower 
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TABLE 1. Total Off-road vehicle trail length as measured from aerial photography maps 
over time in Big Cypress National Preserve. 


Total length Scale of Person/group 
of mapped source who measured 
Year ORV trails map trail length 
1940 160 km 1:63,000 Duever et al., 1979 
1953 250 km 1:63,000 Duever et al., 1979 
1973 1,100 km 1:75,000 Duever et al., 1979 
1984 1,995 km 1:63,000 NPS staff 
1998 46,700 km 1:10,000* Welch and Madden, 1998 


* work is currently being performed to field check result. 


altitudes. Welch and Madden (1998) completed a comprehensive inventory 
of ORV trails that involved magnification of 1:40,000 scale aerial photog- 
raphy taken in 1993 to 1:10,000 scale and mapped all discernable ORV 
features, for a total of 46,700 km of trails. Welch and Madden have more 
recently initiated a field program to verify the accuracy of their estimate. 
The increase in measured trail length over time is summarized in Table 1. 

The trend highlighted by Table | suggests that ORV-induced disturbance 
is occurring faster than the natural system is capable of restoring itself, 
leading to a condition where impacts have accumulated over time. Duever 
and co-workers (1986) noted that once soils are disturbed, few natural mech- 
anisms are capable of restoring the marl soils of the preserve. An analog to 
ORV-induced impacts can be made with past farming operations. Farming 
operations were extensive along the Tamiami Trail and Loop Road prior to 
establishment of the preserve, and occurred mostly in marl soils. There was 
no clearing before cultivation other than burning and furrows were estab- 
lished by modification of the surficial soil (as opposed to rock tilling). De- 
cades later these furrows are still clearly visible (from the air and on the 
ground), suggesting that soil impacts resulting from ORV use will have 
similar persistence over time. 


Wetter conditions of the 1990s relative to the previous two decades— 
Analysis of water level data shows that generally wetter hydrologic condi- 
tions have prevailed in the 1990s relative to the first two decades of the 
preserve’s existence (1970s and 1980s). The degree of disparity between the 
1990s and the previous two decades is illustrated with stage duration curves 
calculated for daily water levels collected at three locations within and ad- 
jacent to the preserve—P34, Bridge 105, and Bridge 84 (Fig. 7). The vertical 
axis reports water level in meters above sea level and the horizontal axis 
reports the time frequency that water levels were at or below the reported 
value. The approximate land-surface elevation is graphically reported for 
each site at the horizontal interface between dark and light background shad- 
ing. Elevation of the land surface does not factor in the spatial variability 
of the surrounding landscape, but is thought to provide a consistent com- 
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Fic. 7. Stage duration curves shown for ten-year intervals of 1970-1979, 1980-1989, 
and 1990-1999 at three long-term monitoring stations in the preserve. Land surface elevation 
is also shown for the purpose of comparing frequency of surface water flooding among the 
three decades. Daily average stage values were used as the baseline data source. 


parison for measuring the duration of surface water inundation among the 
three decades. In the 1990s, standing water is shown to be present 78, 63, 
and 36 percent of the time at P34, Bridge 105, and Bridge 84, respectively, 
in comparison to 25, 30, and 10 percent of the time in the 1970s. Table 2 
reports the duration of surface water inundation at the three stations over 
the past three decades. 

The wetter conditions in the 1990s can be attributed to three factors. (1) 
Increased rainfall: Three years of the 1990s (1994, 1995, 1999) experienced 
rainfall amounts of greater than 50 cm of the average in comparison to no 
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TABLE 2. Summary of surface water flooding duration during the 1970s, 1980s, and 1990s 
at three long-term monitoring stations in Big Cypress National Preserve. Daily average stage 
values were used as the baseline data source. 


Duration of surface water flooding 
(percent of time) 


Station 1970s 1980s 1990s 
Bridge 84 10 8 34 
Bridge 105 30 30 61 
P34 21 71 75 


similar events in the 1970s and 1980s according to the Miami Airport rain- 
fall station (Fig. 8). Drought rainfall seasons of less than 20 cm of the yearly 
average were recorded 10 times in the 1970s and 1980s, with no occurrence 
in the 1990s. (2) Modifications to upstream water management operations: 
Drainage patterns in south Florida were drastically altered in the 1970s by 
the cumulative impacts of physical (canals, levees) and operational (flow 
structures, pump-enhanced drainage, back-pumping) modifications imple- 
mented by the U.S. Army Corps of Engineers (USACE) and South Florida 
Water Management District (SFWMD) from 1947 to the mid 1970s. Mount- 
ing evidence that the new system put too little water into Everglades Na- 
tional Park led to the development of new operational rules and structures 
in the 1980s and 1990s. The wetter conditions documented at P34 in the 
1980s and 1990s are likely a result of modifications to increase flow to Shark 
River Slough in the central park region and to prevent excessive flooding 
of Cape Sable Seaside Sparrow habitat along the east and west banks of the 
slough. (3) Removal of physical barriers or artificial conduits to flow: In- 
creased stage levels recorded at Bridge 84 in the 1990s can be largely at- 
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Fic. 8. Historic rainfall record at Miami International Airport, reported relative to devi- 
ation from the historic average annual year. This graph was modified from a South Florida 
Water Management District poster. 
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tributed to infilling/plugging of the Turner River Canal—a 32 km long canal 
constructed in the late 1950s that was draining interior portions of the 
swamp. 

The persistence of high water in the 1990s is believed to have created 
a condition where soils were inundated with surface water deeper and longer, 
thereby increasing their susceptibility to rutting, trail widening, and estab- 
lishment of new trails from wheeled ORVs. Airboat country has probably 
been less affected, but has created a condition where new trails have been 
developed that were historically too shallow for airboat use. The relative 
dryness of decades preceding the 1990s also leads one to speculate that the 
rise of recreational ORV-use coincided during an era when water levels were 
artificially lower than current levels, and probably caused less visual impacts 
because vast areas were dried down during the hunting season. 

Major restoration activities underway in south Florida, and expected to 
continue for the next several decades, are anticipated to further improve 
natural flow regimes in the central everglades and adjacent ecosystems, in- 
cluding the Big Cypress Swamp. These activities will likely lead to increased 
surface-water stage duration and flow in the eastern portion of the preserve. 


Federal statutes—While only limited guidance was provided in the pre- 
serve’s Enabling Legislation (PL 93-400), the Legislative History of this 
law provides a greater level of detail regarding ORV management. As ex- 
cerpted below, there is an expectation that the preserve will be returned to 
a wilderness character and that ORV use should be managed in tune with 
this principle. 

‘Since the area included in the preserve is largely undeveloped at the 
present time and because it will be managed in a manner which will assure 
its return to the true wilderness character which once prevailed, it will offer 
many outdoor recreation opportunities to the visiting public. During signif- 
icant portions of the year, primitive camping, hiking, and sightseeing will 
be popular activities. Naturalists and wildlife observers will come to see the 
flora and fauna in their natural setting and hunters and fishermen will con- 
tinue to find the area popular. While the use of all-terrain vehicles must be 
carefully regulated by the Secretary to protect the natural, wildlife and wil- 
derness values of the preserve, the bill does not prohibit their use along 
designated road and trails (U.S. Congress, 1974).”’ 

The Legislative History also provided guidance on how preserves, as a 
new unit of the national park system, should be managed in general—by 
merit of Big Thicket National Preserve, Texas, and Big Cypress National 
Preserve, Florida, being the first areas designed as this new unit type in 
1974. The Legislative History states: ‘‘(t)he principal thrust of these areas 
should be the preservation of the natural values which they contain.”’ It also 
states: “‘(n)ational preserves may accommodate significant recreational uses 
without impairing the natural values, but such public use and enjoyment 
would be limited to activities where, or periods when, such human visitation 
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would not interfere with or disrupt the values which the area was created to 
preserve (U.S. Congress, 1974).”’ 

Further guidance regarding ORV management on federal lands is also 
provided in Executive Order 11644, “‘Use of Off-Road Vehicles on the Pub- 
lic Lands,’’ EO 11990. The order states that trails and areas permitting mo- 
torized ORV use must be administratively designated in a way that mini- 
mizes disturbance to natural resources. It also allows for closures or vehicle 
restrictions where adverse impacts are occurring, or where important re- 
sources are threatened. The order also calls for federal agencies to monitor 
the effects of ORVs on lands under their jurisdiction. 


Technological advances since 1974—Technological developments since 
the inception of the preserve have raised new concerns regarding ORV use 
in the preserve. Swamp buggies and airboats have become more powerful 
and more reliable over the past quarter century. In addition, a new type of 
ORV called all-terrain vehicles (ATVs) now offer a low-maintenance, high- 
speed alternative for accessing the backcounty. Cell phones and global po- 
sitioning systems have also become more commonly available. All of these 
developments improve the navigational abilities of ORV operators, allowing 
them to penetrate into unknown territory or previously inaccessible interior 
lands of the preserve. 


CONCLUSIONS—Looking back on 26 years of ORV use in the preserve 
suggests that impacts have been increasing over time. While debate may 
continue on the extent that visual disturbance translates into long-term eco- 
logical damage, the discussion of photographic imagery, trail mapping, hy- 
drologic conditions, legislation, and the advent of new technologies provided 
in this paper are believed to provide strong rationale for modifying the man- - 
ner in which ORVs have been traditionally used in the preserve over the 
past 26 years. Failure to make such changes would cause a continued ac- 
cumulation of disturbance over time, in new areas, and with significant un- 
certainty regarding long-term ecological consequences. Moreover, failure to 
act would be in violation of the laws that were enacted by the United States 
Congress to establish the preserve. 

Big Cypress National Preserve completed a new Off-road Vehicle Man- 
agement Plan in Fall 2000 and is currently in the process of implementing 
it. The plan is believed to meet the challenge of protecting the fragile natural 
resources of the preserve in the spirit of fostering its return to a more wil- 
derness character, while still providing access. The major piece to the plan 
will be the establishment of a designated trail system, ending the era of 
unrestricted ORV access throughout the preserve (Fig. 9). Other components 
of the plan include enacting temporal and spatial closures (heavily disturbed 
areas, Cape Sable Seaside Sparrow habitat, pristine areas), a new three- 
pronged permit system (vehicle permit, ORV operator permit, backcountry 
permit), and initiation of monitoring, research, and restoration. 
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Fic. 9. Map showing off-road vehicle use areas and trails currently being implemented 
in the preserve’s new Off-Road Vehicle Management Plan completed in 2000. Note that the 
illustrated trail network is conceptual. 
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ABSTRACT: Population dynamics and home range in populations of the Key Largo woodrat 
(Neotoma floridana smalli), Key Largo cotton mouse (Peromyscus gossypinus allapaticola) and 
the introduced black rat (Rattus rattus) were investigated in four tropical hardwood hammocks 
on Key Largo, Florida. The hammocks represented three different age classes: medium (dis- 
turbed between 1940 and 1959), medium-old (disturbed between 1926 and 1940), and old 
(disturbed before 1926). N. floridana in medium-old hammock had the highest population size 
and high female reproduction, where R. rattus was absent. N. floridana, in old-aged hammock, 
had low recruitment, low female reproduction, seasonal male reproduction, and a larger home 
range compared to other sites. In the presence of R. rattus, N. floridana had low population 
size, reduced survivorship, and a male biased population. Medium-aged hammocks had the 
highest abundance of P. gossypinus. However, the presence of R. rattus caused P. gossypinus 
to reduce their home range size in those hammocks. R. rattus was only present in a medium- 
aged grid, where it had high population size and recruitment. 


Key Words: Population ecology, Key Largo, Neotoma, Peromyscus, Rattus 


ENDANGERED species can be managed successfully only if their popu- 
lation dynamics and habitat requirements are well understood. Unfortunately, 
the necessary data are not available for many endangered species. For those 
species where data are available, they often consist of nothing more than 
abundance or density in a particular area. Thus, natural resource managers 
may be forced to make critical management decisions on the basis of infor- 
mation. 

Population density may be a poor indicator of habitat quality in small 
mammals (Van Horne, 1983; Ostfeld, 1992). To assess the suitability of a 
habitat, one needs to consider a variety of demographic variables such as 
density, survivorship, recruitment, age structure, reproduction, and space 
use. 

In addition to demography, space use has been shown to vary with 
habitat quality. Home range and movement decreased in food supplemen- 
tation studies of Peromyscus maniculatus (Smith, 1971; Taitt, 1981), Cleth- 
rionomys rufocanus (Ims, 1987), and Microtus townsendii (Taitt and 
Krebs,1981). Animals living in a patchy habitat also should adjust their 
home range and movement accordingly. Diffendorfer and co-workers (1995) 
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reported that movement of three small mammal species increased as habitat 
patchiness increased. ; 

In the tropical hardwood hammocks of Key Largo, Monroe County, 
Florida, high density has been observed to be indicative of habitat quality 
for small mammals by Barbour and Humphrey (1982) and Goodyear (1985). 
The endangered Key Largo woodrat (Neotoma floridana smalli) and the Key 
Largo cotton mouse (Peromyscus gossypinus allapaticola) reside in the trop- 
ical hammock forest of Key Largo that has been characterized as medium- 
to-high productivity rockland hammock (Ross et al., 1992). Along with the 
woodrat and cotton mouse, the black rat (Rattus rattus) occurs in some of 
the hardwood hammocks on Key Largo (Hersh, 1981). 

The greatest threat to both of these endangered subspecies is habitat 
loss. The species were formerly distributed throughout the island but are 
currently found only on the northern half of Key Largo. There are also 
potential threats from exotic species such as black rats, feral cats, and fire 
ants. 

There are approximately 945 hectares (ha) of forest on the northern half 
of Key Largo (Ross et al., 1995). However, not all of this habitat may be 
suitable for woodrats and cotton mice. The greatest densities of woodrats 
and cotton mice have been found in mature hammocks (Barbour and Hum- 
phrey, 1982; Goodyear, 1985). 

Mature hammocks are characterized by little undergrowth, a high degree 
of shade, continuous canopy, and a substantial leaf litter and soil layer. In 
mature hammocks, woodrats nest in the vicinity of tree roots, but will also 
make stick nests if tree root sites are not available. Cotton mice also nest 
around tree roots, in the periphery of woodrat stick nests, and in burrows 
in the leaf litter (Goodyear, 1985). 

Immature hammocks have large amounts of undergrowth with less de- 
veloped root systems, a discontinuous canopy, and a shallow leaf litter. Im- 
mature hammocks are less suitable habitat than mature hammocks for both 
species because of fewer suitable nest sites (Goodyear, 1985; Humphrey and 
Barbour, 1979). 

Hersh (1981) estimated woodrat densities on Key Largo to be 2.2 wood- 
rats per ha in a tract of mature hammock, with little seasonal variation 
between August and May. Humphrey (1988) found 0.3 to 16.4 (mean 7.6) 
woodrats per ha and 3.0 to 44.1 (mean 21.2) cotton mice per ha at six sites 
on Key Largo during the dry season (November through April). He calcu- 
lated totals of 6,500 woodrats and 18,000 cotton mice for the 851 ha of 
habitat he estimated as available. 

All previous studies on the density of the Key Largo woodrat and cotton 
mouse suggest that mature hammock is optimal habitat for both species. 
However, other demographic variables such as survivorship, reproductive 
activity, age structure, and sex ratio were not considered. 

Additionally, Ross and co-workers (1995) classified the hammocks of 
Key Largo into five age classes according to the time since last disturbance. 
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Their five levels of hammock age are: old (disturbed prior to 1926), medium- 
old (disturbed 1926-1940), medium (disturbed 1940-1959), medium young 
(disturbed 1959-1971), and young (disturbed post 1971). Hammock age- 
classes over the entire island was represented as follows: old 27%, medium- 
old 25%, medium 40%, medium-young 2%, and young 6%. This more de- 
tailed classification of the hammocks, in conjunction with population de- 
mography, can now be used to more accurately describe the habitat utilized 
by the small mammal community. 

The aim of this study was to determine population dynamics and space 
use of N. floridana and P. gossypinus on Key Largo in hammocks of dif- 
ferent age classes and compare them to previous studies. Our results on the 
population ecology of these two species should lead to better management 
plans for these two species by providing managers with information on how 
populations of the two native rodents have responded to the protection of 
lands on north Key Largo. 


MeETHODS—This study was conducted on grids in two hammocks on the Crocodile Lake 
National Wildlife Refuge (CLNWR) and two hammocks in the Key Largo Hammocks State 
Botanical Site (KLSBS) on Key Largo, FL from July 1996 to April 1998. Grid 1 was in a 
medium-aged hammock located in CLNWR that covers 16 ha, grid 2 was in a 49 ha old 
hammock in CLNWR, grid 3 was in a 68 ha medium aged hammock in KLSBS, and grid 4 
was in a 29 ha medium-old hammock in KLSBS. Trapping was suspended on grid | in July 
1997 due to an influx of fire ants at the location. 

Medium-young and young age class hammocks were not included in this study because 
they only accounted for only 8% of total hammock area and occurred in patches too small for 
adequate sampling. 

Grids (1.8 ha) were a 10X10 block design with stations spaced 15 m apart. Each station 
had a Sherman XLF-15 (4” X 4.75” X 15’’) live trap. Grids were trapped for three nights 
every four weeks. Captured animals were marked with a uniquely numbered ear tag. On first 
and subsequent captures trap location, body mass, and reproductive condition were recorded. 
Males were considered reproductively active if testes were descended into the scrotum (scrotal). 
Two of the three following characteristics were considered indications of reproductive activity 
for females: perforate vagina, open pubic symphysis, and medium to large nipples size. Age 
class of animals (juvenile or adult) was defined by body mass at which they first became 
reproductive as follows: cotton mice > 22 g, woodrats > 190 g, and black rats > 110 g. 

We used the Jolly-Seber method (Jolly, 1965) to estimate demographic parameters for 
these species. Population parameters estimated for each site were abundance, survivorship, 
recruitment (number of newly tagged animals), sex ratio (% male) of all captures, percent 
reproductive for both sexes, and age structure (% adult). 

Home range size was estimated using the minimum convex polygon method (Mohr, 1947) 
with the program CALHOME (Kia et al., 1996). Ninety-five percent minimum convex polygons 
were used as an estimate of home range size. Only adult individuals with a minimum of five 
captures in successive trapping sessions were included in home range analyses. 

One-way ANOVAs with Tukey’s pairwise comparisons were used to test for differences 
among grids for all variables. Differences of demographic variables for within a grid between 
the wet (May—October) and dry (November—April) seasons were examined using two-sample 
t-tests. All proportions and percentages were arcsine square root transformed (Sokal and Rohlf, 
1995). Data not normally distributed were log transformed. 


RESULTS—Population size—Woodrats and cotton mice were captured on 
all four grids, but black rats only were captured on grid 1. N. floridana 
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Fic. 1. Population size estimates for N. floridana on four grids in tropical hammocks of 
Key Largo, Florida from July 1996 to April 1998. Estimates are Jolly-Seber estimates of pop- 
ulation size. 


abundance (Fig. 1) was significantly different among grids (F = 12.90, d.f. 
= 3, 68, P<0.001) with no differences within grids between seasons (grid 
l:t = 0.76, d.f. = 6, P = 0.48; grid 2:t = —0.22, df. = 16, P = 0.83; grid 
3:t = 0.31, df. = 12, P = 0.76; grid 4:1 = —0.40, d.f. = 10, P = 0.70). 
Tukey’s pairwise comparisons indicated that the population on grid 1 was 
significantly smaller than grids 2, 3, and 4, grid 2 was significantly smaller 
than grid 4, and grid 3 was significantly smaller than grid 4. Mean number 
of individuals + 1 SE on the grids were as follows: grid 1, 1.7+0.3; grid 
2,.3.02250:3= cindy 355247 20s andverid: 465,520) 8: 

P. gossypinus abundance (Fig. 2) also differed significantly among grids 
(F = 10.73, d.f. = 3, 68, P<0.001), with Tukey’s pairwise comparisons 
revealing that grid | had higher populations than grids 2 and 4, grid 2 had 
a smaller population than grid 3, and grid 3 had a larger population than 
grid 4. Populations generally peaked in early winter. Mean number/grid + 
1 SE was as follows: grid 1, 21.1+2.2; grid 2, 14.0+0.9; grid 3, 20.725: 
and grid 4, 14.8 = 0.7. A significant difference in numbers between seasons 
within a grid was found on grid 2 (t = —2.40, d.f. = 17, P = 0.028) and 
grid 4 (t = 3.28, df. = 12, P = 0.006). Grid 2 had higher mean numbers 
in the dry season (15.87 vs. 12.04) while grid 4 had higher numbers in the 
wet season (13.02 vs. 16.86). 

Mean population abundance of R. rattus was 22.9+7.8 individuals. 
Numbers peaked in late fall and declined during the summer, although the 
difference between seasons was not significant. 


Survivorship—N. floridana survivorship was not significantly different 
among grids (Table 1) or within any grid between seasons (grid 1: t = —0.14, 
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Fic. 2. Population size estimates for P. gossypinus on four grids in tropical hammocks 
of Key Largo, Florida from July 1996 to April 1998. Estimates are Jolly-Seber estimates of 
population size. 


eee— 7, F —'0.90; erid-2:¢ = 0.50, df. = 16, P = 0.63; grid 3: t= —0.08, 
df. = 10, P = 0.94; and grid 4: t = —1.6l1, df. = 14, P = 0.13). P. 
gossypinus survivorship was not significantly different between grids (Table 
1) or between seasons (grid 1: t = —1.39, df. = 6, P = 0.21; grid 2: t = 
0.58, df. = 14, P = 0.57; grid 3: t = —0.25, d.f. = 13, P = 0.81; and grid 
4: t = 0.46, d.f. = 16, P = 0.65). Monthly survivorship for R. rattus was 
0.933+0.275 with no significant difference between seasons (t = 0.24, d.f. 
= 3, P = 0.83). 


Recruitment—N. floridana recruitment was not significantly different 
among grids (Table 1) or within any grid between seasons (grid 1: t = —0.27, 
wer — or -— "0.80; srid 2:7 = —0.13, df. = 17; P = 0.90: erid 3: 1 = 
—0.41, df. = 15, P = 0.69; and grid 4: t = —0.62, d.f. = 10, P = 0.55). 
P. gossypinus recruitment was significantly different by site (Table 1) with 
grids | and 2 having significantly greater recruitment than grids 3 and 4. 
There was no significant seasonal difference within any grid (grid 1: t = 
~via. — 7, P — 0.69; end 2: t = —0.50, df. = 13, P = 0.63; grid 3: 
aes, et — 16, P — 0:90; and grid 4: ¢ = —1.25, d.f..= 15, P = 0.23). 
R. rattus recruitment on grid 1 was consistently high with a mean of 
10.7=5.74 new individuals per trapping period, with no significant seasonal 
differences (¢ = —0.41, df. = 6, P = 0.69). 


Age structure—N. floridana age-structure was not significantly different 
among grids (Table 1) or within any grid between seasons (grid 1: t = 0.38, 
oe — 9, P— 071; erid 2: t = —0.46, df. = 20, P = 0.65; grid 3: t = 
—0.20, d.f. = 19, P = 0.84; and grid 4: t = 0.09, d.f. = 16, P = 0.93). In 
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P. gossypinus, grids | and 2 had a significantly lower proportion of adults 
than grid 4. Age structure for this species varied by season on grid 2 (f = 
412, df. = 17, P<0.001) with a mean of 91.6+3.4% adults in the wet 
season and 64.8+5.6% adults in the dry season and grid 4 (t = 3.69, d.f. 
= 16, P = 0.002) with a mean of 96.8+1.7% in the wet season and 
83.6% +3.2 adults in the dry season. R. rattus also had an adult-biased pop- 
ulation (overall mean = 72.1%+11.0) with no difference between seasons 
(¢ = 0.71, d.f. = 10, P = 0.50). 


Sex ratio—N. floridana sex ratio was significantly different among grids 
(Table 1), with grid 1 with a male biased ratio differing from grid 4 with 
female biased sex ratio. There were no significant differences between sea- 
sons within any grid (grid 1: t = 2.40, df. = 6, P = 0.053; grid 2: t = 
oat — 19. P — 0.11; erid 3: t — 0.39, df. = 19, P = 0.70; and grid 
4: t = 0.88, d.f. = 19, P = 0.39). It should be noted that only two females 
(one recaptured once) were caught on grid | during the entire study. 

P. gossypinus sex ratio was also significantly different among grids (Ta- 
ble 1), with grid 4 having a higher proportion of males than grids | and 2. 
Grid 2 was the only grid to have a significant difference in sex ratio between 
seasons (t = 2.39, df. = 20, P = 0.027) with males accounting for 
58.5+3.3% of the captures in the wet season and 46.6 +3.7% in the dry 
season. 

The overall sex ratio of R. rattus on grid 1 was 57.1+8.7% male, with 
no significant difference between seasons (t = —1.08, d.f. = 10, P = 0.30). 


Reproduction—The overall percentage of N. floridana with scrotal testes 
did not differ significantly among grids(Table 1). Only on grid 2 was there 
a significant difference in the percentage of reproductively active males (t¢ 
= 2.08, d.f. = 15, P = 0.05), with a mean of 45.5+16.0% of reproductively 
active individuals in the wet season versus 8.3+8.3% in the dry season. The 
overall percentages of reproductively active N. floridana females differed 
significantly among the three grids for which samples were available, with 
the greatest level of reproductive activity on grid 4 (Table 1). There was no 
difference in reproductive activity between seasons within any grid (grid 2: 
t = 0.00, d.f. = 20, P = 1.0; grid 3: t = —0.08, d.f. = 19, P = 0.94; and 
grid 4: t = 0.63, d.f. = 18, P = 0.53). 

The percentage of reproductively active P. gossypinus males was not 
significantly different among grids (Table 1). Reproductively active males 
occurred throughout the year with a drop in frequency during the late winter/ 
early spring. The seasonal difference was significant on grid 2 with a mean 
of $85.4+£3.3% in the wet season and 63.1+7.1% (@¢ = 2.84, d.f. = 15, P = 
0.012) in the dry season and on grid 3 with a mean of 73.9+2.5% in the 
wet season and 55.8+5.1% (t = 3.17, d.f. = 15, P = 0.006) in the dry 
season. P. gossypinus female reproductive activity did not differ significantly 
among grids (Table 1) or within grids between seasons (grid 1: ¢ = —1.11, 
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TABLE 2. Estimates of home-range size for N. floridana on four grids and overall in 
tropical hammocks of Key Largo, Florida. All data are means + | SE. 


Home range area 


Grid Sex N (m7?) 
Grid 1 Male 1 1575 
Female 0 — 

Grid 2 Male 2, 6413 + 900 
Female ] 1013 

Grid 3 Male 2 213822 450 

Female Z 619 + 169 

Grid 4 Male 5 247 = Os 

Female 6 1463 + 346 

Combined Male 10 3105 + 628 

Female 9 1225-25 257) 


d.f. = 8, P = 0.30; grid 2: t = 0.29, d.f. = 15, P = 0.78; grid 3: ¢ = —0.95, 
d.f. = 15, P = 0.36; and grid 4: ¢ = 0.35, d.f. = 19, P = 0.73). 

For male R. rattus on grid 1, 69.4+10.2% were reproductively active 
with no significant difference between seasons (t = —0.68, d.f. = 9, P = 
0.52). Among females, 57.7+12.3% were reproductively active with no sig- 
nificant difference between seasons (t = 0.90, d.f. = 10, P = 0.39). 


Home range—Male N. floridana had a significantly larger home range 
size than females for all grids combined (F = 7.07, d.f. = 1,17, P = 0.017; 
Table 2). Male ranges also differed significantly among grids (F = 10.34, 
d.f. = 2,6, P = 0.011), with a significantly larger mean size on grid 2 than 
grids 3 and 4. Female home range did not differ significantly among grids 
(F = 1.75, d.f. = 1,6, P = 0.234; Table 2). 

P. gossypinus male home-range was significantly larger (F = 18.55, d.f. 
= 1,106, P<0.001) than that of females (Table 3). Male range size was not 
significantly different among grids (F = 0.76, d.f. = 3,53, P = 0.524), 


TABLE 3. Estimates of home-range size for P. gossypinus on four grids and overall in 
tropical hammocks of Key Largo, Florida. All data are means + 1 SE. 


Home range area 


Grid Sex N (m7) 
Grid 1 Male 8 2672 + 667 
Female ED 947 + 221 
Grid 2 Male 13 3920 + 710 
Female 13 2259 =e 
Grid 3 Male 17 4553 + 876 
Female 18 2288 + 399 
Grid 4 Male 19 3932 + 634 
Female 10 2633 + 590 
Overall Male Sy) 3938 + 383 
Female ais) 2042 + 232 
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although there was a suggestion of smaller home range on grid 1 compared 
with grids 2,3, and 4 (Table 3). Female home-range size (Table 3) differed 
significantly among grids (F = 4.72, d.F. = 3,49, P = 0.006), with grid 1 
having significantly smaller home range size than grids 2, 3, and 4. 

R. rattus on grid 1 had a mean home-range of 1538+582 m°’. R. rattus 
were only captured on grid 1, therefore comparisons among grids were not 
possible. 


DISCUSSION—Previous studies on Key Largo have suggested that the 
most mature hammocks are optimal habitat because of the greater abundance 
of woodrats and cotton mice found there (Barbour and Humphrey, 1982; 
Goodyear, 1985). Hersh (1981) estimated woodrat densities on Key Largo 
to be 2.2 woodrats per ha in a tract of mature hammock, with little seasonal 
variation between August and May. Humphrey (1988) found 0.3 to 16.4 
(mean 7.6) woodrats per ha and 3.0 to 44.1 (mean 21.2) cotton mice per ha 
at six sites on Key Largo during the dry season (November through April). 
Humphrey (1988) sampled in a variety of hammocks of different ages, but 
did not describe those ages. He noted a north-south gradient with woodrats 
and cotton mice being less abundant in the northern part of Key Largo. 

Densities in our study ranged from 0.9 to 3.1 (mean 1.8) woodrats per 
ha and 7.8 to 12.1 (mean 9.9) cotton mice per ha with no north-south gra- 
dient as reported by Humphrey (1988). The results presented here for wood- 
rats are consistent with those reported by Hersh (1981). These differences 
may be due to natural fluctuations in population size over time or reflective 
of the quality of the habitat sampled. 

In our study, N. floridana population size was greatest in medium-old 
hammock, intermediate in medium and old hammock, and lowest in a me- 
dium hammock where R. rattus was present. Recruitment was highest in 
medium-old hammock, which also had greater female reproduction with a 
female-biased sex ratio. The medium hammock with R. rattus had the lowest 
survivorship, a low proportion of adults, and a population that was essen- 
tially all male. 

The old and medium hammock without R. rattus had high survivorship 
and intermediate population density, with female reproductive activity higher 
in the medium than in the old hammock. The higher percentage of females 
in reproductive condition suggests that medium hammock is of higher qual- 
ity than old. Also, home-range in the old hammock was greater than on 
other age classes which also suggests that old hammock is lower quality 
habitat. 

For P. gossypinus population size was highest in medium aged ham- 
mock. Although both medium age hammocks had high population size, re- 
cruitment was more variable in the hammock where R. rattus was present. 
The high abundance of individuals in medium hammocks suggests that this 
age class is optimal habitat for P. gossypinus. The low density and female 
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reproductive activity on grid 2 suggests that old hammock was the sub- 
optimal for P. gossypinus. 

The population dynamics of N. floridana and P. gossypinus on grid 1 
were potentially affected by the presence of R. rattus. The proximity to 
human development of grid | may account for the abundance of R. rattus. 
It should be noted that R. rattus was reproductively active individuals were 
recaptured in the hammock and are not a floating population of non-repro- 
ductive animals as reported by Humphrey (1992). N. floridana had the low- 
est survivorship, lowest proportion of adults, and a highly male-biased sex 
ratio on grid | (Table 1). Although both N. floridana and P. gossypinus were 
possibly affected by the presence of R. rattus, the extremely poor perfor- 
mance of N. floridana is of particular concern as R. rattus has been impli- 
cated in the extinction of several rodents throughout the Greater and Lesser 
Antilles (Woods, 1989). The results presented here suggest that P. gossy- 
pinus can coexist with R. rattus far better then N. floridana. 

If intermediate post-disturbance stages of tropical hammock are the best 
quality habitat for P. gossypinus and N. floridana, as suggested by the data, 
this indicates that periodic hurricanes may play a role in maintaining opti- 
mum habitat quality for these species. Duever and co-workers (1994) stated 
that the Florida Keys experience a major hurricane every 27 years. Distur- 
bance by hurricanes will result in the necessary heterogeneous matrix of 
habitat allowing source populations of both species to persist and colonize 
suitable hammock stages during post-hurricane recovery. 

Two other factors potentially affecting populations of these species on 
Key Largo are the presence of R. rattus and fire ants (Solenopsis wagnert). 
The latter appear to be colonizing new areas on Key Largo. Fire ants were 
not present on grid | at the start of the study, but became so abundant that 
trapping could not continue on the area. 
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USING GEOGRAPHIC INFORMATION SYSTEMS TO 
EVALUATE AVAILABLE HABITAT FOR TWO RODENT 
SPECIES ON NORTH KEY LARGO, FLORIDA. 


BRITT KEITH AND MICHAEL S. GAINES 


University of Miami, Department of Biology, Cox Science Building Room 215, 
1301 Memorial Drive, Coral Gables, FL 33146 


ABSTRACT: We used Geographic Information Systems (GIS) as an assessment tool to 
predict the amount and distribution of available habitat for the Neotoma floridana smalli (Key 
Largo woodrat) and the Peromyscus gossypinus allapaticola (Key Largo cotton mouse) on North 
Key Largo, Florida. Three times as many cotton mice were live trapped in the hammock than 
were woodrats. None of the habitat variables used in the analyses; hardwood hammock age 
class, fragment size, runway density and presence of exotic vegetation, were significantly as- 
sociated with Key Largo woodrat abundance estimates. Key Largo cotton mouse abundance 
was highest in medium-old hammocks, smaller fragments, hammocks with high runway densities 
and in hammocks where the exotic vegetation was absent. A stepwise multiple regression anal- 
ysis revealed that hammock age class and fragment size explained most of the variation in 
cotton mice abundance. ArcView’s Spatial Analyst GIS was used to analyze the stepwise mul- 
tiple regression equations of the four habitat variables. The output was a map of the habitat 
which predicted areas of high (greater than 30) and low (less than 30) cotton mice abundance. 


Key Words: GIS, habitats, Neotoma floridana smalli, Peromyscus gossy- 
pinus allapaticola 


IN 1994 the Florida Game and Fresh Water Fish Commission and the 
United States Fish and Wildlife Service listed Neotoma floridana smalli 
(Key Largo woodrat) and Peromyscus gossypinus allapaticola (Key Largo 
cotton mouse) as endangered. These two small mammal subspecies have 
been of concern to wildlife biologists because they are in danger of becom- 
ing extinct due to habitat destruction and possible competition by invading 
exotic species. 

The Key Largo woodrat and the Key Largo cotton mouse are the south- 
ernmost subspecies of their respective genera (Brown, 1970) and are com- 
pletely isolated on North Key Largo, Florida (Fig. 1). Habitat requirements 
of the two species in North Key Largo are sketchy. Among others, a major 
factor that may be influencing the habitat requirements of these animals is 
the maturity level of the hardwood hammocks. Brown (1978 a, b) claimed 
that woodrats and cotton mice are limited to mature hardwood hammocks, 
and that both species are absent in medium to young age stands. However, 
Goodyear (1985) has trapped both species in both mature and immature 
hammocks. 
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Enlarged area 


_ Hardwood hammock 


Fic. 1. Map showing the study site of North Key Largo. 
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The objective of this study is to utilize Geographic Information Systems 
(GIS) to identify where the most suitable habitat occurs for woodrats and 
cotton mice in North Key Largo. 


METHODS AND MATERIALS—Study area—All fieldwork took place within the tropical hard- 
wood hammocks in North Key Largo. There are no continuous stretches of undisturbed ham- 
mock, due to roads and development in the area, therefore the study area was divided into nine 
fragments that varied in size from 25-85 ha. 


Trapping methods—We determined the distribution of woodrats and cotton mice by con- 
ducting systematic line transect sampling over the entire extent of the hardwood hammock in 
North Key Largo. The total number of transects for the study site was 35. Each transect was 
sampled once. Stations were flagged at 10-m intervals along a transect for a total of 25 stations 
per transect. Two ventilated Sherman XLF-15 small mammal live traps were placed at each 
station. The traps were opened and baited in the evenings, and checked for occupancy in the 
early mornings. Each animal caught was ear tagged for future identification using individually 
numbered Monel ear tags. 


Abundance estimates—The Schnabel method of abundance estimation (Krebs, 1989) was 
calculated for both species at each transect. Ross and co-workers (1995), classified the hardwood 
hammock into five different age classes related to when they were last disturbed. These five 
classes are old (undisturbed since 1926), medium-old (disturbed 1926-1940), medium (dis- 
turbed 1940-1959), medium-young (disturbed 1959-1971) and young (disturbed since 1971). 
These age classes were the primary habitat features used for the analysis. Because the habitat 
is severely fragmented, a regression analysis was performed using fragment size as the inde- 
pendent variable and the abundances of each species as the dependent variable. 


Runway densities—Anecdotal evidence during field trapping indicate that animals use 
runways for escape. Runways are defined as fallen branches, logs, vines, rocks, or any com- 
bination of the four. The density of runways was recorded in the field at every line transect. 
An analysis of variance (ANOVA) was used to detect any statistically significant differences 
of runway densities on the transects in the five hammock age classes. In addition, a regression 
analysis was performed with runway densities as the independent variable and the abundances 
of each species as the dependent variable. A regression analysis was also performed with 
fragment size as the independent variable and runway density as the independent variable. 

Goodyear (1985) reported that exotic vegetation might have a negative in this area. There- 
fore, we incorporated exotic vegetation as a structural feature of the habitat. For the analysis, 
we used presence or absence, without considering the densities of exotics. 


GIS methods—Arc View GIS is the software used throughout this study. All transects were 
mapped using a Trimble GPS unit. This transect layer was overlaid onto an existing digital line 
map of North Key Largo. In addition, layers and their corresponding attribute tables were 
created for each habitat variable used in this study. Hammock age class, fragment size, runway 
density and the presence or absence of exotic vegetation were potential candidates for GIS map 
construction. A stepwise multiple regression analysis was performed using these candidates as 
the independent variables, and the abundance data as the dependent variable. The resulting 
regression equations were input directly into ArcView to create a suitability map of the habitat, 
ranked according to the predicted abundance of animals in different areas. 


RESULTS—Woodrats were captured on 21 of the 35 transects. ANOVA’s 
and regressions found no statistically significant differences in woodrat 
abundance among hammock age classes (F = 1.59, df = 3,30, p = 0.211), 
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Cotton Mouse Abundance 


1 2 3 


Hammock Age Class 


1 = old hammock 
2 = medium-old hammock 
3 = medium hammock 


Fic. 2. Bar graph showing the relationship between cotton mouse abundance and ham- 
mock age class. 


fragment size (F = 0.24, df = 1,5, p = 0.64. R* = 0.046), runway density 
G@ — 0.10, df — 1,5, p = 0.760, R? = 0.020) or whether or not exotic 
vegetation was present on or near the transects (F = 3.28, df = 1,33, p = 
0.079). 

Cotton mice were captured on 27 of the 35 transects. Three times as 
many cotton mice were captured than woodrats (153 as opposed to 47). An 
ANOVA detected a statistically significant difference of cotton mice when 
hammock age classes were lumped, favoring higher abundance estimates in 
medium-old hammock over old and medium (F = 6.39, df = 2,5, p = 0.042, 
Fig. 2). There was a statistically significant negative relationship when cotton 
mice abundance was regressed as the dependent variable on fragment size, 
the independent variable (F = 6.25, df = 1,7, p = 0.041, R* = 0.472). There 
was also a Statistically significant difference favoring higher cotton mice 
abundance estimates where exotic vegetation was absent on the transects (F 
= 4.87, df = 1,33, p = 0.034). 

Runway densities varied across the site and between each hammock age 
class. When cotton mice abundance was regressed on runway density there 
was a Statistically significant relationship (F = 5.33, df = 1,7, p = 0.054, 
R? = 0.432). There were statistically higher runway densities occurring in 
the medium-old hammock age class, compared to medium, old and young 
hammocks respectively (F = 3.93, df = 3,29, p = 0.018). There was also 
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a Statistically significant negative relationship when fragment size was re- 
gressed on runway densities (F = 7.14, df = 1,7, p = 0.032, R? = 0.505). 

Unfortunately, a habitat suitability map for the woodrats could not be 
created because of small sample sizes. The ANOVA’s for the different in- 
dependent habitat variables (hammock age class, fragment size, runway den- 
sity and presence/absence of exotics) showed no significant differences in 
woodrat abundance estimates. However, significant results for the cotton 
mice provided sufficient information for a GIS map using the habitat vari- 
ables. An estimated 790 ha of the total 945 ha was classified as suitable for 
the cotton mice. A stepwise multiple regression using all habitat variables 
approached statistical significance (F = 5.22, df = 4, 4, p = 0.069, R* = 
0.678). This regression equation was used when creating the GIS model due 
to its potential biological significance. The stepwise multiple regression re- 
vealed that hammock age class and fragment size were the only statistically 
significant variables in the model. Together they explained 76.6% of the 
variation (F = 9.83, df = 2, 6, p = 0.013). Figure 3 shows the results of 
the stepwise multiple regression GIS models, showing predicted areas of 
high (>30) and low (<30) cotton mice abundance. 


DIscUSSsION—There have been only a few studies documenting the hab- 
itat requirements of the Key Largo woodrat and the Key Largo cotton mouse. 
The literature states that woodrats are more specific in their habitat require- 
ments than cotton mice (Goodyear, 1985; Brown, 1970). Our results indicate 
the reverse. There was no clear association of habitat variables with woodrat 
abundance. This could have been due to their small sample size and that the 
habitat variables that we chose may not be ecologically relevant to the or- 
ganisms. However, cotton mice did show significant differences in abun- 
dance for different habitat variables. Thus, we will focus on cotton mice. 

Several authors have claimed that cotton mice prefer the most mature 
hammocks. Our results do not support claims in the literature that cotton 
mice prefer only the oldest hammocks. In our study, cotton mice were 
trapped on 77% of the transects. Medium-old, old and medium hammock 
age classes all had a statistically significantly higher abundance of cotton 
mice than medium-young and young age classes. Runway densities are also 
significantly higher in medium and medium-old age classes compared to old 
hammocks. Higher runway densities were also positively correlated with 
higher abundances of cotton mice. One possible explanation for this is that 
the medium age class probably provides enough cover for these animals, 
but the canopy allows enough light penetration for the growth of additional 
plants and vines which provide runways. Old hammocks have a closed can- 
opy, therefore an open understory and less variability. 

In addition to the individual regressions with the habitat variables, the 
stepwise multiple regressions revealed some interesting results. When all of 
the variables were included, the model approached statistical significance. It 
is the co-linearity of these variables that may complicate which variables 
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Fic. 3. A map of the hardwood hammock in North Key Largo, produced by the stepwise 
multiple regression analysis, predicting areas of hammock which have a low (O—30) and high 
(30-60) abundance of cotton mice. 
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most affect the animals. Although hammock age class explains the most 
variation in cotton mouse abundance (56%), the issue is whether these an- 
imals prefer the oldest hammocks, which is assumed in the literature. This 
study reveals that the medium old age classes had the highest cotton mice 
abundances. 

The final GIS map highlighted areas in the hardwood hammock where 
the models predict different cotton mouse abundances. There are approxi- 
mately 945 ha of hammock on North Key Largo. The map identified 790 
of the 945 total ha of hammock deemed suitable for cotton mice. These are 
areas where the suitability of the hammock increases as the values of cotton 
mice abundance also increase. Thus, a large percentage of the hardwood 
hammock on North Key Largo is deemed suitable at varying degrees for 
cotton mice, according to the GIS model (Fig.3). It takes into account all of 
the habitat variables that were initially predicted as important in a statistical 
sense to these animals. If this model were to be validated in the field, the 
areas represented in black should yield an abundance of at least 30 cotton 
mice per transect. This model was tested against the abundance estimates 
found in this study, yielding a 68% concordance with the data. That is, 68% 
of the areas predicted to have at least 30 cotton mice in fact did, according 
to the original trapping data collected. 

Geographic Information Systems are quickly evolving to produce more 
finite and accurate results, particularly in the field of ecology. At present, 
the technology works best on large-scale landscapes where many habitat 
types may be classified, including a plethora of variables. Although this 
North Key Largo study was done on a smaller scale, it still provides useful 
information to agencies trying to manage this parcel of land. We are hopeful 
the map generated from this study will help to identify areas of special 
concern. This may allow the agencies to more easily customize the guidance 
given to the biologists working to protect and upgrade the ecosystem and 
endangered species in North Key Largo. 
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ABSTRACT: The alligator snapping turtle, Macroclemys temminckii , is capable of detect- 
ing and responding to the presence of chemosensory cues from two naturally—occurring prey 
species, the Mississippi mud turtle, Kinosternon subrubrum, and the common musk turtle, S. 
odoratus. Preference tests were conducted in an experimental chamber in which juveniles of 
M. temminckii were exposed to dechlorinated water in one-half of the chamber and an odor 
plume from water that had been occupied by either S. odoratus or K. subrubrum. Juveniles of 
M. temminckii spent significantly more time in water containing prey odor cues and also 
exhibited more rapid movements and more frequent episodes of gular pumping, as compared 
to those exhibited by turtles in untreated water. The ability to utilize odor cues is discussed 
with respect to detection of prey and patch residence time. 


Key Words: Chemosensory cues, Macroclemys, predator recognition, turtles 


PREVIOUS research has shown that many animals can reduce the risk of 
predation by utilizing chemosensory cues associated with potential predators 
(see review by Kats and Dill, 1998). Most of the research has focused on 
certain invertebrate species (Chivers et al., 1996) as well as fish (Hara, 
1992), larval amphibians (Petranka et al., 1987; Manteifel, 1995) and mam- 
mals (Kats and Dill, 1998; Barreto and Macdonald, 1999). Among the rep- 
tiles, the use of odor cues to detect predators has been reported for lizards 
(Cooper, 1990; Phillips and Alberts, 1992; Van Damme and Castilla, 1996) 
and snakes (Burger, 1989). Turtles, in contrast, have received far less atten- 
tion (Jackson, 1990). 

Less well known, but equally important, is the ability of predators to 
utilize chemosensory cues to locate potential prey. This has been well es- 
tablished in many species of mammals (Galef et al., 1994) and carrion- 
feeding birds (Kare and Mason, 1986). Punzo and Kakoyi (1997) showed 
that wolf and lynx spiders can learn to recognize odor cues associated with 
naturally-occurring prey species, and furthermore, the presence of these cues 
increases patch residence time for these spiders. 
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Although most species of turtles possess a well-developed visual capac- 
ity (Mrosovsky and Boycott, 1966), they are also capable of utilizing che- 
mosensory cues to locate preferred microhabitats (Yeomans, 1995; Quinn 
and Graves, 1998) and food (Burghardt and Hess, 1966; Hatanaka and Mat- 
suzaki, 1993). Jackson (1990) showed that musk turtles (Sternotherus minor 
peltifer and S. carinatus) utilized chemical cues to detect and avoid a nat- 
urally-occurring predator, the alligator snapping turtle (Macroclemys tem- 
minckii) . The purpose of this study was to determine whether or not M. 
temminckii can utilize chemical cues to detect the presence of musk and 
mud turtles. 


MEeETHODS—AII turtles used in this study were obtained from eggs hatched in captivity. 
The eggs were deposited by females originally collected at several locations in Livingston 
Parish, LA, in July, 1999. Eggs were obtained from 5 females of M. temminckii (clutch size 
ranged from 12-31 eggs), 8 females of Kinosternon subrubrum (Mississippi mud turtle; 2—5 
eggs per female), and 6 females of Sternotherus odoratus (common musk turtle; 2-8 eggs per 
female). Alligator snapping turtles are known to feed on musk and mud turtles (Allen and Neill, 
1950; Jackson, 1990). 

After collection, eggs were incubated in a Havabator incubator (Bush Herpetological Sup- 
ply, Neodosha, KS) at 27°C and 72% relative humidity. After hatching, juvenile turtles were 
maintained individually in 7-gallon aquaria containing dechlorinated water, artificial plants, a 
gravel substrate, and a rock on which the animals could climb to get out of the water. Water 
in the aquaria was changed twice per week. Turtles were maintained on a diet of Reed’s Turtle 
Food (Bush Herpetological Supply) and mealworms (Tenebrio molitor) . All turtles were 2 
weeks old when tested and they were all naive (i.e.—none had prior exposure to any other 
turtle prior to testing). 

Experimental trials were conducted in a stainless steel chamber (135 80 X< 35 cm) based 
on the design described by Jackson (1990). To summarize, the chamber was filled with 105L 
of dechlorinated water replenished at a rate of 1OL / min. A 95 X 25 cm piece of clear plastic 
was affixed to the floor of the chamber, standing on its long edge at the upstream end at the 
midline of the tank width. This effectively divided the tank into 39.4 cm wide halves, except 
at the downstream end of the chamber. The inflow of water to each half of the chamber entered 
through | mm holes in a capped plastic pipe and was directed into a submerged trough which 
was buffered with plastic dowels so that the entering water did not cause any turbulence. Water 
exited through 6 mm diam. holes in the downstream end wall. The addition of colored food 
dye verified the presence of laminar flow characterized by a distinct boundary where the flows 
of each half were in contact downstream from the plastic divider. In several test runs, no color 
could be detected in the chamber 10 min after adding the dye. All observations were recorded 
with a Panasonic MV5O0 video camera from behind a one-way mirror to minimize disturbance 
to the animals. 

At the beginning of each trial, an individual M. temminckii was given a choice between 
plain dechlorinated water (control) or dechlorinated water containing 2% of water from an 
aquarium that had housed 5 individuals of either S. odoratus or S. subrubrum for a 3-day period 
(treatments). Before water reached the experimental chamber each perforated trough also re- 
ceived 100 ml / min of water via gravity flow from either the control or one of the two 
treatments. The flows of treatment and control water were stopped at the end of each trial, but 
mass flow of water continued for at least 20 min between trials. Control and treatment inputs 
were renewed 20 min before the next trial was initaiated. 

Each trial started by placing one M. temminckii, facing upstream, at the midline of the 
chamber halfway between the downstream end of the chamber and the downstream end of the 
plastic divider. The movements of the turtle were recorded and quantified for a 20-min period 
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TABLE 1. The time (as proportion of total), mean speed (cm/min), and frequency of epi- 
sodes of gular pumping (EGP) of Macroclemys temminckii in the two halves of the experimental 
chamber as a function of the type of stimulus (dechlorinated water vs. prey species: Kinosternon 
subrubrum and Sternotherus odoratus). Data expressed as means (N = 45). Values in paren- 
theses represent (+SD). See text for details. 


S. odoratus K. subrubrum 
Parameter Water Prey Water Prey 
Time 0.14 (+0.06) 0.86 (+0.23) OF 2071) 0.88 (+0.26) 
Speed 15.3 (#4.2) DIGEST) IS Ges y7)) 32.4 (+4.9) 
EGP 12 Oriiy) Df KEZ2D) 0.9 (+0.04) 6:2) (2231) 


using a grid coordinate system (comprised of l-cm square grids) affixed to the walls of the 
chamber. In this way, we were able to determine the amount of time spent in treated vs. 
untreated halves of the chamber and the speed at which the turtle moved (cm / min). We also 
recorded the presence or absence of gular pumping by each subject since this behavior is known 
to increase as a function of locomotor activity (Hatanaka and Matsuzaki, 1993). 

Two replicates were conducted on 25 and 20 M. temminckii during these experiments. 
Subjects were chosen randomly from all of the hatchling turtles. We utilized a completely 
randomized block design. Each turtle received one trial with one of the two treated water 
supplies as an alternative. The proportions of time spent in each half of the chamber were 
subjected to an arcsine transformation and did not deviate significantly from normality (Sokol 
and Rohlf, 1995). When variances were almost equal, two-tailed t tests were used to determine 
differences between means. When variances were not equal, we used Welch’s approximation 
to the t test (Sokal and Rohlf, 1995). We measured speed of movement by calculating the 
number of cm a turtle moved per min. Data from turtles that were inactive or were absent from 
either half of the chamber were not used in the determination of speed. 


RESULTS AND DIsCUSSION—AII of the juvenile alligator snapping turtles 
exhibited ambulatory movements within the experimental chamber. The re- 
sults shown in Table 1 indicate that M. temminckii did not spend equal 
amounts of time in treated vs. untreated halves of the chamber; rather, this 
predator spent significantly more time in the area of the chamber containing 
cues from the prey species (P < 0.004 for both S. odoratus and S. subru- 
brum). Individuals of M. temminckii also moved more quickly in the pres- 
ence of prey chemosensory cues (P < 0.01 for both prey species). 

The frequency of episodes of gular pumping by M. temminckii also 
increased significantly in the presence of prey cues (Table 1) as shown by 
a two-tailed t test (P < 0.001 for both prey species). In many trials (N = 
34), individual alligator snapping turtles in the area of the chamber contain- 
ing dechlorinated water would begin to move toward the boundary of the 
water containing prey cues and immediately increase their pumping of the 
gular region before proceeding into the treated water. 

These experiments demonstrate that alligator snapping turtles are capa- 
ble of detecting and responding to the chemosensory cues of musk and mud 
turtles. In its natural habitat, M. temminckii is a cryptic predator that typically 
chooses a patch and remains motionless waiting for prey to pass close by 
(Allen and Neill, 1950). However, it also moves slowly along the bottom of 
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ponds, lakes, and rivers searching for food (Ashton and Ashton, 1985). It 
uses its worm-mimicking tongue to lure fish, but can also sieze turtles and 
crustaceans that come within striking distance (Ernst and Barbour, 1972). 
Based on the results of this study, M. temminckii may utilize chemosensory 
cues from potential prey species to make decisions concerning patch choice 
and patch residence time. An earlier study by Burghardt and Hess (1966) 
on the common snapping turtle, Chelydra serpentina, showed that this spe- 
cies is capable of responding to the odors of certain foods that it experiences 
early in life, resulting in a preference for those food types later in life. 
Similar findings have been reported for some insects (Kats and Dill, 1998), 
spiders (Punzo and Kakoyi, 1997), and other vertebrates (Burger, 1989; Coo- 
per, 1990; Hara, 1992; Galef et al., 1994). 

The increased rate of gular pumping exhibited by M. temminckii in the 
presence of prey odor cues is in agreement with previous research suggesting 
that aquatic turtles utilize gular pumping to enhance their ability to chemi- 
cally sample the surrounding water (Manton, 1979). 

In conclusion, predators utilize chemosensory cues to detect the presence 
of prey species or habitats (patches) recently occupied by prey (Punzo and 
Kukoyi, 1997; Cooper, 2000). Likewise, prey species utilize chemosensory 
cues to avoid potential predators (Manteifel, 1995; Chivers et al., 1996; Kats 
and Dill, 1998). Musk turtles are capable of detecting the presence of M. 
temminckii by utilizing chemosensory cues associated with this predator 
(Jackson, 1990). This study shows that M. temminckii can, in a similar fash- 
ion, detect the presence of musk and mud turtles. 

The detection of chemosensory cues can also influence the rate of move- 
ment of prey and predators. Prey animals frequently stop moving or move 
at slower speeds after they have encountered predator odor cues (Jackson, 
1990; Weldon, 1991). In this study, predators increased their rate of move- 
ment (appetitive behavior) in the presence of odor cues associated with prey. 

Because these experiments were conducted on naive juvenile turtles, the 
preference for water containing musk turtle chemosensory cues may be in- 
nate in this species. Future studies should attempt to determine whether this 
preference can be modified based on past experience as well, and whether 
or not M. temminckii can detect and respond to odor cues from other prey 
species. 
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ABSTRACT: Typically total lead has been taken as a measure of the toxicity of lead in a 
sample, whereas a direct measurement of toxicity would seem to be a more meaningful measure. 
This study summarizes analysis of lead using a Microtox 500 toxicity analyzer. The instrument 
examines the effect of samples on the luminescent bacterium Vibrio fischeri. We compared the 
effect of standard solutions of aqueous lead nitrate with soil samples containing known amounts 
of lead nitrate. The LD5;, for soluble lead was about 2 ppm. There was no notable difference 
between a control sample (lead-free soil) and high-organic soil containing about 1000 ppm 
lead. The LD, for both samples was about 50,000 ppm. The latter value is likely a function of 
the nature of the soil type and the organic content, but it does demonstrate that high-organic 
soil is capable of absorbing soluble lead. It also demonstrates the need for a specific mea- 
surement of toxicity as a guide to bioavailability and bioaccumulation. 


Key Words: Bioaccumulation, humic acids, lead, Microtox, soil, toxicity 


LEAD compounds, their speciation, and their distribution, constitute a 
widespread health problem in urban areas (Hafen and Brinkmann, 1996; 
Mushak, 1992, U.S. Senate, 1990; Mielke et al., 1989). Lead-based disorders 
are thought to be the most common environmentally caused health problems 
in the United Sates (CDC, 1985; U.S. Senate, 1990). Samples of lead com- 
pounds are found in such common sources as lead-based paints, lead-based 
solder in water pipes, and automobile exhaust particulates in soils (Hafen 
and Brinkmann, 1996). A study of lead in soils adjacent to an interstate 
highway near the USF Tampa campus found that over one-third of the sam- 
ples contained lead concentrations (>500 wg/g) considered to be hazardous 
by the Environmental Protection Agency (Hafen and Brinkman, 1996). 

One aspect that deserves more consideration is the speciation of lead 
and the bioavailability of the lead in a given sample. In short, if the lead 
were present as the lead equivalent of apatite, Pb,(PO,),OH, a highly insol- 
uble form of lead, it is unlikely that its bioavailability would be very high. 
Comparatively little attention has been given to organically bound lead, de- 
spite the greater toxicity (Grandjean and Grandjean, 1984), and such species 
have been found in wine (Lobinski et al., 1994). And thus, the toxicity of 
lead samples may be of as much interest as the total lead concentration. 

One of the key problems of toxicity testing often has been the need for 
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a reliable, rapid method. One such method is the use of a Microtox analyzer, 
which exposes a luminescent bacterium, Vibrio fischeri, to aqueous samples 
and measures the change of light output by the test organism (cf. Hao et al., 
1995, 1996; Guzella et al., 1996; Chen et al., 1997; Devalls et al., 1997). 
The bacterium is grown under optimal conditions by the same firm that 
produces the instrument. When needed, the lyophylized bacterium is rehy- 
drated with a reconstitution solution to provide a suspension of organisms. 
The analyzer measures the light output of the bacteria after they have been 
exposed to a sample and compares this number to the light output of a 
control (reagent blank). The difference in the light output between the test 
sample and the control is attributed to the effect on the organism. Each 
cuvette contains approximately a million test organisms, and the system 
measures the simultaneous light output of all the test organisms. 

This study describes the use of a Microtox analyzer to assay the toxicity 
of lead in water and soil samples and the results that were obtained. 


MATERIALS AND METHODS—Reagents and materials—Lead nitrate (analytical grade) was 
obtained from J.T. Baker. Nitric acid (69.6%), hydrogen peroxide (30%), and hydrochloric acid 
(35%) were reagent grade chemicals purchased from Fisher Scientific. Soil was commercial 
potting soil from K-Mart and was stored in Pyrex containers. 

Samples of Vibrio fischeri were obtained from Azur Environmental and maintained in 
freeze-dried form at —30°C until ready for use, where upon a given ampoule was treated with 
Azur Reconstituting Solution, as directed. Diluent was added to the sample and control cuvettes, 
and Osmotic Adjusting Solution was added to the control cuvette. 


Instrumentation—A Microtox 500 analyzer, obtained from Azur Environmental, was cou- 
pled to a Gateway 2000 computer and Hewlett Packard 610C printer. The Microtox instrument 
measures both acute and chronic toxicity; in these studies, acute toxicity was measured, using 
a 5-min study period. 


Analysis of soil samples—Water content was determined as loss on drying at 70 °C for 
24 h to constant weight. The weight loss was taken as the weight of moisture and the percentage 
weight loss is 100 times the weight loss divided by the weight of the initial sample. Organic 
content was determined by the Ball procedure (Ball, 1964). The dried sample in a tared con- 
tainer was placed in a muffle furnace, and the temperature raised to and held at 375°C for 17 
hours. Organic carbon was taken as the loss in weight of a dried sample (Eqn. 1). 


Percent organic carbon = (loss in weight/weight of original sample) < 100 (1) 


Lead content—Soil was mixed with an aqueous solution of lead nitrate, the mixture was 
stirred for four hours, and dried in an oven (70°C for 24 hours). Lead content was determined 
using atomic absorption spectrometry, following an extraction procedure (McGrath and Cun- 
liffe, 1985; Hafen and Brinkmann, 1996). Specifically, a 2-g sample of soil in a 150 mL beaker 
was treated with 10 mL of dilute nitric acid (8 M) and heated at 95°C for 15 min. The solution 
was cooled, and 5 mL of conc. nitric acid (16 M) was added and heated for an additional 30 
minutes at 95°C. The last step was repeated and the volume was reduced to about 5 mL by 
evaporation. The sample was cooled, and treated with 2 mL of deionized water and 3 ml of 
30% hydrogen peroxide. The sample was heated gently to initiate the reaction, then 5 mL of 
conc. HCl (16 M) and 10 mL of deionized water were added, and the sample warmed to near 
boiling. The sample was then allowed to cool to room temperature and was filtered through 
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TABLE 1. Effect of lead nitrate in water in five minutes, using a Microtox 500. 


Lead %Effect* on the bacteria Vibrio fischeri (Smin-test) 
concen- 
tration Dilution Dilution Dilution Dilution 
(ppm) (5.625%) GuiE25 2) (22.50%) (45%) 
10 17.47 29.95 46.26 68.12 
9 20.74 41.17 53.68 VAM 
8 9.362 31.86 48.39 72.24 
7 1339 26.75 47.09 64.53 
6 S295) 8087 46.03 Horo 
5 ESS 14.75 34.75 59.62 
4 14.87 DNS 54.54 
3 —** — 13.64 Bel 
2 — 10.07 2523.0 41.41 
1 Scaae = ae; = 


* J bacteria killed. 
** __: No valid results were obtained. 


Whatman No. 42 filter paper. Samples were then analyzed using a Varian atomic absorption 
spectrometer, equipped with a lead lamp. 


Toxiciry TESTING—General procedure—An Azur Microtox 500 unit was to test the tox- 
icity, and the so-called “‘wizard’”’ test or guideline procedure (part of the instrumental software) 
was used to decide which Microtox test was suitable for a given sample type. Acute toxicity 
of lead was measured by a basic test in these experiments. 

The test sample was aqueous lead nitrate, which was prepared from an amount of reagent 
grade lead nitrate in deionized water to obtain 100 ppm lead standard solution (0.000483 M). 
Then the standard solution was diluted to series solutions. 


RESULTS AND DISCUSSION—Data_ treatment—Following the Microtox 
procedures gives percentage effect as a function of percentage dilution for 
a test of a given duration (typically 5 minutes). Take for example aqueous 
lead nitrate using a five-minute test (Table 1). The table lists four serial 
dilutions (45, 22,5, 11.25, and 5.625 %) of a given concentration (ranging 
from 1—10 ppm Pb), and summarizes percentage of the organisms that were 
destroyed (for a given dilution of a particular initial concentration of aqueous 
lead nitrate). In addition, a longer reaction time (15 min.) was also used to 
test for toxicity of lower concentrations (Table 2). 

Next, probit analysis (Finney, 1971) was used to calculate the effective 
concentration needed to kill 50% of the sample, LD). This process trans- 
forms a S-shape survival curve to a linear plot, and the LD. is easier to 
obtain from the linear plot because a probit of 5 corresponds to a 50% kill 
or survival. The percentage killed values were converted to probits, using 
an available table (Finney, 1971), and the results are given (Fig. 1). Indi- 
vidual concentrations could be tested, but pooling the data (N = 34) gave 
an LD,, of 1.9 ppm, which was consistent with the values obtained for 
smaller numbers of samples of given concentrations. 
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TABLE 2. Effect of lead nitrate in water in fifteen minutes, using a Microtox 500. 


%Effect* on the bacteria Vibrio fischeri (15min-test) 


Lead con- 
centration Dilution Dilution Dilution Dilution 
(ppm) (5.625%) (11.25%) (22.50%) (45%) 
10 48.17 70.67 83.16 90.47 
9 83.87 89.93 90.72 90.89 
8 76.53 88.59 90.63 92.02 
7 72.97 85.89 90.40 91.36 
6 71.48 84.85 89.86 90.12 
5 60.97 79.41 89.09 91.6 
4 Die JMS S 99.96 90.33 
3 rk Se 75.24 84.52 
2 6.821 d027, 80.68 88.87 
1 — — 9.036 11.76 
* % bacteria killed. 
** __: No valid results were obtained. 


Lead analyses—homogeneity—lIt is necessary to make sure that the soil 
sample was homogeneous with respect to lead distribution before analyzing 
for toxicity. Accordingly, ten samples were extracted with 1M nitric acid 
solution, a more rapid, but less efficient, method than the EPA-based method. 
Using the simpler procedure, the mean lead content was found to be 5.16 


=34 


Probit 
(ee) 


3 y = 1.5471x + 4.5658 
R? = 0.9122 


-0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 
Log concentration 


Fic. 1. Probit analysis of all Microtox data (N = 34) for effect of aqueous lead solution 
on V. fischeri. Here, y = 1.5471 + 4.5658; r? = 0.9122; when Y = 5 (50% kill), LD50 = 1.9 
ppm Pb*=. 
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TABLE 3. Toxicity of original soil for five- and fifteen-minute runs using a Microtox 500. 


5-min % Effect* on 15-min % Effect on 

Sample Soil conc., ppm Vibrio fischeri Vibrio fischeri 
Control 0 

1 386.7 10.89 0.485 

2 773.4 11.82 1.805 

eS) 1547 337) 3.800 

4 3094 Ike 9.760 

5 6188 8}.5)1| 9.655 

6 12380 29.79 24.17 

a 24750 3)9)-3))! SAT) 

8 49500 55.80 56.96 

9 99000 74.21 59.02 


* JY bacteria killed. 


ppm with a relative standard deviation of 5.1%. This indicated the homo- 
geneity of the sample. Unfortunately, the calculated value was 6.24 ppm, so 
the mean percent extraction was 82.7%. Subsequently, using the EPA-based 
method, two analyses gave 90 and 97.8% extraction. 


Toxicity of lead samples—The data for toxicity of lead-containing soil 
were treated similarly to the data for aqueous lead. The basic solid-phase 
test was used. This allows contact between the bacterial medium and the 
solid sample and provides quantitative data on solid samples. 

First, untreated soil samples were tested because of the assumption that 
soil could have toxic substances present. Second, treated samples of soil 
(containing 996 ppm Pb) were studied and the effect was measured as a 
function of time (Table 4). The results appear to be consistent for 5-, 15-, 
and 30-minute exposures and fairly consistent results ;were obtained for the 
greater concentrations and the results were similar to those obtained for plain 


TABLE 4. Toxicity of lead-treated soil (996 ppm Pb) for five- and fifteen-minute runs using 
a Microtox 500. 


5-min % Effect*  15-min % Effect 30-min % Effect 


Sample Soil conc., ppm on Vibrio fischeri on Vibrio fischeri on Vibrio fischeri 

Control 0) 
1 386.7 6.080 3.040 1.380 
2, VIBE: 2.040 33505 9.740 
3 1547 7.400 9.040 1.590 
4 3094 9.170 16.25 5.850 
a 6188 nips ial 37'S) 
6 12380 2220 29.20 31.10 
7 24750 38.60 44.65 45.56 
8 49500 44.95 56.08 55.48 
9 99000 SO) MS) 61.95 68.12 


* Y pacteria killed. 
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(untreated) soil. Compare, for example, the data for 49,500 ppm soil for 
which the results at 5- and 15 min tests were 56% and 57% respectively, as 
compared with the results for lead treated soil for which results were 46% 
(S-min), 56% (15 min) and 55% (30 min). 

The most significant observation concerning the effect of soil on the 
test bacteria show the results for the two different soils (without and with 
lead) is that there is only slight difference between the two. This conclusion 
supports the view that this amount of soluble lead (996 ppm) in this high- 
organic soil does not differ in toxicity from untreated soil. In contrast, 2 
ppm soluble lead in water was the LD, for the test bacteria. 

It should not be greatly surprising that soluble lead would become un- 
available in soil with a high organic content. Other metals could have been 
present; the point made here is that the soluble lead became unavailable, 
and there was no difference between test and control for treated lead. High- 
organic soils are known to have significant amounts of chelating enties that 
can entrap metal ions (Langford et al., 1983; Tuschall and Brezonik, 1983; 
Weber, 1983). What would be surprising would be to discover such high- 
organic soils very widely distributed in Florida, which seems to be covered 
with a fine-grained quartz sand of low mineral content. 
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FROM THE EDITORS—Those Who Make It Possible—II 


When you look at this copy of your journal, it may seem obvious who 
is responsible—The authors and, oh yes, the Editors. True enough, but the 
truth doesn’t go far enough in this instance. Reviewers, dedicated volunteers, 
willing to share their expertise and time are vital, as we try to acknowledge 
annually. In addition, when you see the mention of Allen Press below the 
masthead (on the inside front cover), most readers don’t realize the large 
number of persons there who help us and strive for a quality product. These 
include Ms. Suzanne Wood and her colleagues at Allen Press. Also on the 
masthead is the name of a major player, Dr. Richard Turner, the dedicated, 
able Business Manager, who worries about money matters and reminds one 
of the revised truism, ‘““We aren’t in business to make money, but we must 
make [recover some] money to remain in business.”’ From time to time, Dr. 
Richard Turner will arrange for advertisements that will appear in an issue 
and help provide the financial support that is a must for continued operation; 
one ad appears in this issue. We hope that more such supporters will come 
forward. 

In supporting roles are other important players—Dr. Gay Biery-Hami- 
liton and Ms. Natalie Smith in the Executive Secretary’s office keeping track 
of members and sending out reminder notices. We have valued the leadership 
of presidents, past and present. Most recently, we have been delighted by 
the support of President Barry Wharton, whose leadership role led to the 
Academy web page and to obtaining space for the Florida Scientist abstracts 
on APT On Line, a web page connected to our web page. This is a web 
site where interested persons can efficiently search Abstracts and authors so 
that work published in the Florida Scientist becomes more accessible and 
more widely used. And it was President Wharton who encouraged the listing 
of some five past years of the Florida Scientist. In this he was supported 
by members of the Academy Council. And finally, since the journal is an 
official publication of the Academy, we are grateful to all the members, a 
portion of whose dues make it possible to continue. And so to all who are 
involved in the publication of Florida Scientist, our sincerest thanks for your 
support. 
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4) sshoppers 
of Florida 


The first volume in the 
Invertebrates of Florida series! 


Grasshoppers of Florida 
John L. Capinera, Clay W. Scherer, 
and Jason M. Squitier 

“Florida naturalists will be delighted 
by Grasshoppers of Florida. This won- 
derfully presented book unlocks an 
amazing world that previously was 
unavailable to most backyard and 
weekend naturalists.”—Richard 
Franz, Florida Museum of Natural 
History 


This first field guide to the grasshop- 
pers of Florida provides identification 
of the 70 species known to occur in 
Florida. The authors present the 
biology, behavior, ecological signifi- 
cance, and damage potential of grass- 
hoppers and emphasize their interre- 
lationship with Florida’s fauna. A 
straightforward, detailed description 
of each species is accompanied by a 
distribution map. Nearly all species 
are pictured in color, and detailed line 
drawings help distinguish among 
those more difficult to identify. 


176 pp. 102 color plates, 42 drawings, 70 
distribution maps. Flexibind, $34.95 


Order through full-service 
booksellers, our website at 
www.upf.com, or with 
VISA or M/C toll free: 
1-800-226-3822 


UNIVERSITY PRESS OF 


The Cuban Treefrog 

in Florida 

Life History of a Successful 
Colonizing Species 

Walter E. Meshaka, Jr. 

“The abundance of data cited in this 
study not only characterizes the 
Cuban Treefrog but establishes it as a 
model of a successful colonizer. This 
latter feature vastly increases the 
importance of this work by providing 
wildlife biologists with a list of char- 
acteristics that can be applied to other 
introduced species to determine the 
likelihood of successful establishment 
and expansion in non-native habi- 
tats—particularly those severely 
altered by human activity.”—Robert 
Powell, Carnegie Museum of Natural 
History, author of A Key to the Amphib- 
ians and Reptiles of the Continental 
United States and Canada 


Walter Meshaka presents the natural 
history of the Cuban Treefrog from the 
perspective of its phenomenal success, 
in terms of sheer numbers and geo- 
graphic range, as a colonizer in South 
Florida and, in particular, the Ever- 
glades. It supplies a methodology for 
evaluating and setting priorities for 
the threats facing Florida’s amphibian 
and reptile populations and identifies 
the most vulnerable species, provid- 
ing a base for management decisions. 
It also presents and interprets a large 
data set associated with patterns of 


colonization and predictions. 
208 pp. 83 figures, 71 tables. Cloth, $69.95 


André Michaux in Florida 


An Eighteenth-Century Botanical 


Journey 

Walter Kingsley Taylor 

and Eliane M. Norman 

“The appeal of André Michaux in 
Florida lies in its interest and useful- 
ness to both professional and amateur 
botanists, horticulturalists, gardeners, 
field naturalists, and historians.” 
—Ronald L. Stuckey, Museum of 
Biological Diversity, Ohio State 
University 


March. 288 pp. 30 b&w photos, 16 maps. 
Cloth, $39.95 
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Life History of a Successful Colonizing Species 


Watter E. Meshaka Jr 


El Nino in History 
Storming Through the Ages 


César N. Caviedes 

“A staggering amount of historical 
research to pinpoint El Nino events in 
time and space, interpret them using 
physical theory, and correlate them 
with human history.”—Betty J. 
Meggers, Smithsonian Institution 


César Caviedes provides the first 
comprehensive historical account of El 
Nino. Combining scientific accuracy 
with readable presentation, he brings 
together all existing information, 
references, and clues about past El 
Nino occurrences and their impact on 
political, military, social, economic, 
and environmental history. 


304 pp. 4 tables, 64 drawings, 12 color photos. 
Cloth, $24.95 


Florida Hurricanes 


and Tropical Storms 
1871-2001, Expanded Edition 


John M. Williams and Iver W. Duedall 
From reviews of the first edition: 

“A comprehensive, factual, concise 
book on Florida hurricanes. .. . Land 
and space photographs and storm 
trajectory maps are especially infor- 
mative.”—Choice 

June. 176 pp. 76 b&w photos, 20 color photos, 


4 tables, 13 tracking charts. 
Paper, $14.95 


Gainesville, Tallahassee, Tampa, Boca Raton, Pensacola, Orlando, Miami, Jacksonville, Fort Myers 
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INSTRUCTION TO AUTHORS 


This information is available at two web sites: 
(1) IES site: http:www.cas.usf.edu (click on “Centers and Institutes”, then select “Institute for Environmen- 
tal Studies”, then select “Florida Scientist’). 


(2) FAS site: www-floridaacademyofsciences.org (select “Florida Scientist’). 


It is possible to search for papers and abstracts of papers for the past six years: Go the FAS site (above), select 
“Allen Press”. When the link is complete, select “Advanced Search” and follow the directions. 


The Florida Academy of Sciences Web page, indicated above, allows interested persons to join the Florida 
Academy of Sciences and pay dues. 
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